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I.  INTRODUCTION 


Theories  to  predict  penetration  of  shaped  charges  have  met  with  varying 
degrees  of  success.  Except  for  relatively  short  standoffs,  penetrations 
diverge  greatly  and  never  really  approach  idealized  values.  Inspection  of 
jets  through  the  use  of  flash  radiographs  reveals  why  idealized  penetration 
is  not  a  reality:  because  the  idealized  jet  is  not  a  reality.  The  presence 
of  drift  velocity  (seen  in  a  radiograph  as  jet  curvature)  is  suggested  as 
being  the  main  cause  for  the  variance  in  jet  penetrations.  Transverse  (or 
alternately  drift)  velocity  refers  to  the  motion  of  shaped-charge  jet  particles 
in  the  plane  perpendicular  to  the  charge  axis  of  the  warhead  (drift  or  trans¬ 
verse  plane).  The  distance  travelled  by  a  jet  particle  in  the  transverse 
plane  is  referred  to  as  the  particle  deviance. 

All  shaped  charge  jets  contain  drift  velocity  to  some  degree  or  another. 
The  origins  of  drift  are  naturally  linked  to  imperfections  in  the  warhead 
itself.  Lack  of  uniformity  of  liner  thickness,  nonhomogeneous  explosive,  as 
well  as  asymmetric  liner  confinements  must  all  be  considered  as  possible 
causes  of  drift  velocity  in  any  given  round.  Fissures  in  the  explosive 
material  and  uncentered  boosters  can  also  contribute  to  jet  performance 
that  is  less  than  perfect.  With  respect  to  the  rounds  studied  in  this 
analysis,  it  was  also  noticed  that  rounds  with  pressed  explosive  possessed 
consistently  lower  penetration  capabilities  than  their  cast  explosive  counter¬ 
parts.  One  theory  is  that  the  high  pressures  (20,000  psi)  required  to  press 
the  explosive  were  capable  of  distorting  the  liner.  The  problems  charac¬ 
terizing  these  pressed  rounds  were  typical  of  the  early  days  of  pressing. 
Technology  in  the  area  of  pressing  explosives  has  since  improved  to  a  point 
where  many  of  the  problems  then  inherent  in  the  rounds  studied  have  been 
corrected.  Some  uncased  pressed  rounds  now  show  penetrations  of  5.5  calibers 
at  24  caliber  standoff.  Nonetheless,  the  pressed  rounds  that  were  studied 
possessed  adverse  qualities  that  manifested  themselves  in  poor  jets.  Thus, 
with  a  multitude  of  sources  for  asymmetric  conditions  always  present,  it  is 
not  surprising  to  find  the  regularity  with  which  drift  velocity  occurs  in  jets. 

An  algebraic  analysis  can  be  performed  to  evaluate  numerically  the  drift 
velocities  of  the  particles  of  a  given  jet.  Preliminary  and  somewhat  simpli¬ 
fied  computations  indicate  that  drift  velocity  alone  is  sufficient  to  account 
for  major  variations  in  jet  penetration. 


II.  MEASURING  THE  DRIFT  VELOCITY  OF  A  PARTICLE 

The  flash  radiograph  is  the  best  tool  to  perform  detailed  analysis  of 
the  shaped  charge  jet  in  flight.  It  would,  therefore,  be  appropriate  to 
explain  the  nature  and  setup  of  the  apparatus  and  how  it  relates  to  the  data 
obtained  and  used  in  the  analysis.  A  diagram  of  the  apparatus  is  shown  in 
Figure  1.  The  jet  is  radiographed  from  three  different  angles  at  three 
different  times  before  it  finally  penetrates  the  target.  The  angular  sepa¬ 
ration  of  each  of  the  x-ray  tubes  is  45®.  A  sample  radiograph  of  one  of 
the  rounds  studied  is  shown  in  Figure  2.  Note  the  fiducial  markings  which 


Figure  1.  Radiographic  Apparatus  Used  to  Record  Flight  of  Shaped-Charge  Jet 


2.  Sample  Radiograph  of  a  Shaped-Charge  Jet  in  Flight 


lie  not  on  but  parallel  to  the  charge  axis.  Data  are  extracted  from  the 
radiographs  by  a  digitizing  process  whereby  each  particle  on  the  film  is 
defined  by  six  points  located  on  the  perimeter  of  the  particle.  By  this 
means,  calculations  are  made  for  each  particle's  length,  diameter,  velocity, 
mass,  momentum,  and  kinetic  energy  using  procedures  described  by  Blische 
and  Simmons . 1 

If  it  is  assumed  that  the  fiducials  are  parallel  to  the  charge  axis  of 
a  round,  the  radiograph  of  an  ideal  jet  would  contain  an  image  of  the  jet 
parallel  to  the  fiducials  on  all  exposures.  This  assumption  is  not  strictly 

true,  but  shall  be  discussed  in  Section  IV.  With  near  perfect  jets  rarely 

being  observed,  it  becomes  important  to  measure  the  transverse  movement  of 
particles  with  respect  to  the  charge  axis.  However,  the  location  of  the 
charge  axis  is,  in  itself,  unknown.  One  assumes  that  the  charge  axis  is 

parallel  to  the  fiducials  on  all  radiograph  exposures,  though  it  is  not 

coincident  with  any  one  of  them.  It  becomes  necessary,  therefore,  to  choose 
a  workable  reference.  The  rear  of  the  jet  has  always  been  considered  the 
most  stable  portion.  For  this  reason,  a  particle  near  the  rear  of  the  jet 
that  is  readily  definable  in  all  three  radiographs  is  chosen  as  the  reference. 
The  reference  axis  is  an  axis  parallel  to  the  warhead  charge  axis  on  which 
the  reference  particle  lies.  If  the  reference  particle  is  not  drifting  in 
the  transverse  plane,  the  reference  axis  will  be  coincident  with  the  charge 
axis . 


To  see  how  the  transverse  velocity  of  a  particle  might  be  measured,  it 
would  be  convenient  to  imagine  looking  down  the  charge  axis,  as  if  from 
liner  to  target.  Ignoring  all  particles  of  the  jet  except  the  reference 
particle  and  the  particle  in  question,  the  charge  axis  view  of  the  situation 
could  be  described  by  Figure  3.  Note  that  the  film  shown  in  the  figure  is 
just  the  "edge"  of  the  film  since  the  view  is  down  the  charge  axis.  The 
deviance  (perpendicular  distance  between  the  reference  axis  and  the  particle 
in  question)  that  is  seen  in  the  radiograph  is  less  than  or  at  most  equal 
to  the  actual  deviance  due  to  the  cosine  effect .  Given  the  apparatus  of 
Figure  1,  and  having  exposed  the  three  radiographs  simultaneously,  the 
reference  particle  and  the  particle  in  question  could  be  precisely  located 
with  respect  to  each  other  in  the  transverse  or  drift  plane.  This  situation 
is  illustrated  in  Figure  4.  The  time  delay  that  actually  exists  among  the 
three  exposures  complicates  matters  little  and  three  governing  equations  are 


t  *  w 

readily  derived: 

*  D  /cos(0+45°) 

VT  — - 

* 

(1) 

1 


v*  D_/cos0 

VT  "  — t— - 

(2) 

t2 

JJ.  Bti&che,  8.  Simon*,  "A  M eXhod  (>o>i  Reducing  Vcrta  uom  Radiograph*  o{, 
Shaped- Charge  Jet*,"  SRL  Report  ARBRL-TR-02330,  Jun  81  (At?  At 02770). 
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Figure  3.  View  of  Particle  Deviation  as  Projected  onto  the  Transverse  Plane 


REFERENCE 

PARTICLE 


Locating  a  Particle  in  the  Trans' 


Dj/cos(0-4S°) 


where  -  apparent  deviance  read  off  the  K'th  radiograph 

t^  -  the  time  between  particle  formation  and  the  K'th  exposure 

0  -  deviation  angle  between  path  of  deviation  and  some  arbitrary 

reference  (chosen  parallel  to  film  2) 

Equations  (1),  (2),  and  (3)  all  contain  two  unknowns,  namely  magnitude 
of  drift  velocity  (V*)  and  angle  of  deviation  (0) .  Any  two  of  these  three 

equations  can  be  solved  simultaneously  to  yield  a  result  for  transverse 
velocity,  both  in  magnitude  and  direction.  In  reality,  the  equations  are 
solved  in  all  three  combinations,  and  a  weighted  average  is  taken  depending 
on  the  particle's  orientation  with  the  radiographs.  This  velocity  shall  be 
called  by  the  name  of  relative  drift  velocity  since  the  particles  are  moving 
with  respect  to  a  moving  reference  (the  reference  particle). 


III.  THE  DRIFT  VELOCITY  DISTRIBUTION  FOR  A  GIVEN  ROUND 

When  relative  transverse  velocities  are  evaluated  for  all  the  particles 
of  a  jet,  a  polar  plot  can  be  made  whereby  a  line  segment  with  length  pro¬ 
portional  to  the  magnitude  of  drift  velocity  is  plotted  at  the  appropriate 
deviation  angle  in  the  transverse  plane.  A  sample  velocity  fan  is  shown  in 
Figure  5.  This  relative  velocity  fan  shows  a  kind  of  spray  of  the  jet  parti¬ 
cles  away  from  the  reference  particle  in  the  deviation  plane. 

When  a  plot  is  made  of  the  axial  versus  transverse  particle  velocity 
for  a  given  round,  a  graph  similar  to  Figure  6  frequently  results.  The  data 
seems  to  be  of  the  form: 

VA  =  ^T  *  VREF 


where:  V^  -  axial  velocity  of  particle 

* 

V.j.  -  relative  transverse  velocity  of  particle 
K  -  constant  varying  from  round  to  round 
V„.  -  axial  velocity  of  the  reference  particle 

Kfcr 

Most  rounds  studied  were  in  high  correlation  with  Equation  (4) . 

The  next  quantity  to  determine  is  the  drift  velocity  of  the  reference 
particle  from  which  all  the  calculations  were  made.  If  we  assume  that  the 
reference  drift  velocity  is  negligible,  then  the  separation  distance  between 
impacts  of  the  jet  tip  and  the  reference  particle  would  routinely  exceed 
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Figure  5.  Relative  Velocity  Fan  Showing  the  Pattern 
Formed  by  the  Drifting  Jet  Particles 
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1/2  charge  diameter  (CD)  at  a  24  CD  standoff  for  typical  drift  distributions 
studied.  Under  such  conditions,  the  reference  particle  .rill  never  enter  the 
crater  created  by  the  jet  tip  because  of  excessive  drift.  However,  it  is 
not  uncommon  for  a  jet  to  produce  a  single  entrance  hole  into  the  target  at 
24  CD  standoff.  Thus,  the  reference  must  possess  a  drift  velocity  which  will 
tend  to  bring  the  reference  into  the  crater  formed  by  the  jet  tip.  An  equation 
is  needed  which  expresses  that  the  reference  particle  should  deviate  the  same 
distance  from  the  charge  axis  as  the  jet  tip  upon  reaching  a  target  block  at 
a  given  standoff.  Mathematically,  the  assumption  is: 


where : 


Dr  -  the  distance  the  reference  has  deviated  upon  reaching  the 
target  surface  (magnitude  only) 

D.  -  the  distance  the  jet  tip  has  deviated  upon  reaching  the 
J  target  surface  (magnitude  only) 


(5) 


Equation  (5)  does  not  preclude  drift  among  the  jet  particles  since  it 
pertains  to  magnitude  only.  An  angular  separation  of  the  drifting  particles 
still  exists  and  will  tend  to  separate  the  jet  particles  at  extended  standoff. 
For  a  small  angular  deviation  among  jet  particles  however.  Equation  (5) 
gives  jets  the  ability  to  produce  single  crater  target  impacts.  It  is  sug¬ 
gested  solely  because  a  compensatory  reference  particle  drift  was  required  and 
is  the  best  which  can  be  proposed  at  this  time.  Since  drift  velocity  is  assumed 
to  be  constant  after  a  particle's  formation.  Equation  (5)  can  be  expressed  as: 

\  lR  ■  VT  <6> 


where:  VT  -  the  absolute  drift  velocity  of  the  associated  particle 

t  -  the  time  required  for  the  associated  particle  to  reach  the 
target  surface 

R  -  the  reference  particle 

j  -  an  arbitrary  jet  particle 

But  the  time  required  for  a  particle  to  reach  the  target  surface  is  equal 
to  the  standoff  as  measured  from  the  virtual  origin  (ZQ)  divided  by  the 

particle's  axial  velocity  (V^J .  Therefore,  Equation  (6)  expands  to: 


Simplifying: 


(7) 


(8) 
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A  RELATIVE  VELOCITY  FAN  WITH 
SCALED  REFERENCE  VELOCITY  ARROW 
INDICATING  MAGNITUDE  AND 
DIRECTION  OF  REFERENCE  DRIFT.  THE 
APEX  OF  THE  FAN  REPRESENTS  THE 
REFERENCE  PARTICLE  FROM  WHICH 
DRIFT  MEASUREMENTS  WERE 
ORIGINALLY  TAKEN. 


AN  ABSOLUTE  VELOCITY  FAN  IS 
CREATED  BY  CONNECTING  THE  TAIL 
OF  THE  REFERENCE  VECTOR  TO  THE 
RELATIVE  VELOCITY  FAN.  THE  APEX 
OF  THE  FAN  LOCATES  THE  CHARGE 
AXIS  IN  THE  TRANSVERSE  PLANE. 
THE  REFERENCE  PARTICLE  NOW 
BECOMES  JUST  ANOTHER  DRIFTING 
PARTICLE. 


Figure  7.  Transfonnation  of  Relative  Fan  Into  Absolute  Fan 
for  a  Hypothetical  Drift  Distribution 


Figure  8.  Effect  of  Transformation  from  Relative  t> 
Drift  Velocity  on  the  Slope  of  the  Curve 


For  a  crater  profile  in  which  all  particles  of  the  jet  enter  the  target  sur¬ 
face  at  the  sane  point.  Equation  (S)  nust  hold  for  all  particles,  implying: 


tj-  =  constant  (9) 

T 

Since  the  reference  is  at  the  apex  of  the  relative  transverse  velocity 
fan,  the  reference  velocity  would  appear  as  a  vector  with  its  tip  attached  to 

the  apex  of  the  relative  velocity  fan.  Its  length  corresponds  to  the  drift 

velocity  of  the  reference  particle  and  its  direction  indicates  the  direction 
in  which  the  reference  is  deviated.  Its  direction  of  deviation  is  assumed  to 

be  about  the  same  as  that  of  the  slowest  digitized  particle.  The  absolute 

velocities  of  the  balance  of  the  particles  are  calculated  by  vectorially 
adding  the  reference  velocity  to  the  relative  velocities  of  the  respective 
particles.  Thus,  an  absolute  drift  velocity  fan  is  created  by  drawing  the 
segments  from  the  tail  of  the  reference  vector  to  the  tips  of  the  relative 
drift  velocity  fan. (See  Figure  7.)  Graphically,  Equation  9  implies  a  linear 
VA  vs.  VT  plot  which  intercepts  the  origin.  With  the  proper  vectorial  addi¬ 
tion  of  the  reference  drift  as  just  described,  relative  drift  distributions 
like  Figure  6  readily  map  into  absolute  drift  distributions  which  satisfy 
Equation  9. (See  Figure  8.) 

There  is  no  way  to  know  whether  this  transformation  into  absolute  veloc¬ 
ities  is  the  correct  one.  However,  it  is  believed  that  the  directional 
assumption  of  the  reference  is  a  safe  one,  and  that  the  Equation  (5)  assump¬ 
tion  is  as  good  as  any  considering  the  need  for  the  drift  correction  of 
particles.  In  addition,  the  results  obtained  from  the  transformation  are 
supported  by  the  experimental  data  •  Note  that.  Equation  (9)  does  not  guarantee 
that  all  particles  will  only  form  a  one -crater  hole  profile  even  though  all 
the  particles  should  deviate  the  same  distance  from  the  charge  axis.  A  large 
angular  spread  in  a  jet  (characterized  by  a  wide  absolute  velocity  fan)  may 
produce  enough  relative  deviation  among  the  particles  so  as  to  cause  multiple 
impacts  at  long  standoff.  Thus,  the  deviating  factors  in  jets  are  both  the 
angular  spread  and  the  magnitudes  of  the  drift  velocities. 


IV.  THE  QUESTION  OF  TILT 

A  possible  and  potentially  fatal  flaw  in  the  previous  discussion  concerns 
the  presence  of  tilt.  The  word  tilt  is  taken  to  mean  a  charge  axis  that  is 
not  parallel  to  the  radiograph  fiducials.  An  original  assumption  was  that 
there  was  no  tilt,  or  that  charge  axis  and  fiducials  were  aligned.  It  is 
obvious  that  some  sort  of  tilt  must  exist,  however  small.  The  answer  as  to 
when  tilt  becomes  critical  is  a  bit  disconcerting.  Depending  on  the  jet,  a 
one-degree  tilt  of  the  charge  axis  off  the  fiducials  could  induce  an  apparent 
drift  velocity  of  160  m/sec  in  the  jet  tip  as  measured  with  respect  to  the 
fiducials.  It  is  likely  that  uncertainties  of  this  magnitude  exist  regularly. 


The  consequence  of  tilt  on  a  jet  is  an  apparent  drift  velocity  distribu¬ 
tion,  which  is  also  governed  by  Equation  (9),  which  is  the  assumed  relation  for 
actual  drift  velocity.  One  may  rightfully  wonder  whether  all  of  the  drift 
velocities  read  off  the  radiographs  are  due  completely  to  tilt,  and  whether 
graphs  like  that  of  Figure  6  are  a  consequence  of  tilt  alone.  In  fact,  what 
is  read  off  the  radiographs  is  the  superposition  of  both  tilt  and  real  drift 
velocity  distributions.  However,  it  should  be  noted  that  an  important  aspect 
of  tilt  is  that  it  cannot  produce  a  velocity  fan.  All  of  the  particles  of  a 
jet  under  the  influence  of  tilt  exclusively  would  appear  to  drift  in  a  uniform 
direction.  Yet,  Equation  (9)  has  been  seen  to  hold  for  jets  with  velocity 
fans  spanning  in  excess  of  90°. 

The  concern  here  is  whether  a  nonlinear  drift  distribution  could  be 
"tilted"  into  appearing  linear  on  the  radiographs,  and  thus  explain  the  data 
without  having  the  actual  drift  obey  Equation  (9).  Such  could  only  be  the 
case  if  the  magnitude  of  actual  drift  were  small  in  relation  to  tilted  drift. 

If  this  were  the  case,  however,  the  angular  spans  of  the  velocity  distributions 
would  also  be  quite  small,  which  is  usually  not  the  case.  Thus,  it  appears 
that  the  actual  drift  velocities  are  at  least  of  the  same  order  of  magnitude 

as  the  apparent  drift  velocities  due  to  tilt.  Thus,  it  seems  that  the  actual 

drift  distributions,  though  probably  different  from  the  measured  drift  dis¬ 
tributions,  can  still  be  adequately  described  by  Equation  (9) .  This  belief 
is  based  only  on  the  rounds  studied.  No  explanation  is  offered  as  to  why  there 

might  be  a  relation  between  the  axial  and  transverse  velocities  of  jet  parti¬ 

cles.  Also,  it  is  uncertain  as  to  whether  this  may  be  a  general  phenomena  or 
characteristic  only  of  the  rounds  studied. 

Support  for  Equation  (9)  was  found  in  work  done  by  Aerojet  General 
Corporation. 2  Aerojet  was  interested  in  the  crater  formation  of  individually 
impacting  particles.  A  shaped  charge  supplied  the  particles  for  their 
experiment.  The  method  employed  for  dispersion  of  the  particles  was  asymmetric 
initiation.  Considering  that  it  was  desired  to  induce  much  drift  velocity  in 
the  particles,  the  effects  of  tilt  could  be  neglected  in  this  experiment.  The 
target  plate  used  in  their  experiment  is  shown  in  Figure  9.  The  particle 
impacts  generally  lie  at  a  constant  distance  from  the  slug  impact.  Assuming 
the  slug  flight  path  to  be  coincident  with  the  charge  axis.  Equation  9  is  the 
governing  equation. 

A  penetration  code  was  created  to  predict  penetrations  for  nonideal 
shaped-charge  jets.  Unlike  the  DSM  model,3  which  can  express  penetration  in 
a  single  equation,  the  model  developed  analyzes  the  jet  piecemeal.  It  is  des¬ 
cribed  in  detail  in  the  following  section. 


?  K.  N.  Kneyenhagen,  J.  E.  Fengason,  R.  R.  Randall,  J.  P.  Joyce,  "Special 
Explosive  Pnojectons Pnoc.  6th  Symposium  hypenveloclty  Impact,  Vol.  I 
[Held  Cleveland,  Ohio,  Apnll  30 -May  2,  7963). 

a 

R.  VIPenslo,  J.  Simon,  A.  Menendlno,  "Penetnatlon  o£  Shaped- Change  J ets 
Into  Metallic  Tangets,"  Ballistic  Reseanch  Labonatony  RepoAt  No.  1296, 
September,  1965  (AP  476717). 
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Figure  9.  Target  Plate  for  Asymmetrically  Initiated 
Shaped  Charge  (Reference  2] 
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V.  THE  CREATION  OF  A  PENETRATION  CODE  TO  ANALYZE  DRIFT 


The  purpose  of  analyzing  drift  velocity  distributions  of  shaped-charge 
jets  is  to  see  how  these  distributions  affect  penetration.  The  dynamics  of 
penetration  is  not,  however,  a  thoroughly  understood  subject.  Many  theories 
exist  in  attempts  to  discribe  the  mechanical  properties  of  matter.  Classical, 
empirical  and  even  wave  mechanical  concepts  are  used  in  pursuit  of  the  subject. 
Consequently,  rather  than  adding  to  the  spectrum  of  theories  that  already  exist, 
the  author  chose  to  draw  upon  information  already  available  on  the  subject 
whenever  possible  and  to  combine  views  in  a  seemingly  compatible  manner. 

One  shaped-charge  penetration  model  that  has  received  much  attention  in 
recent  years  was  formulated  by  DiPersio,  Simon,  and  Merendino  (the  DSM  model).3 
This  model  does  a  good  job  at  predicting  penetration  of  shaped  charges  even  at 
extended  standoffs.  The  model  accomplished  this  however,  by  implying  that 
there  exists  a  cutoff  penetration  velocity,  a  minimum  penetration  velocity  below 
which  a  jet  cannot  penetrate.  Unfortunately,  the  DSM  cutoff  penetration 
velocity  implies  that  jet  material  traveling  below  certain  hypervelocity 
speeds  is  incapable  of  penetrating  target  material.  Substantial  quantities 
of  jet  material  are  thus  eliminated  from  the  penetration  process.  The  DSM 
model  is  empirically  based  and  really  attempts  only  to  describe  the  data  at 
hand  and  does  not  enter  an  intense  theoretical  discussion  to  describe  the 
origins  of  jet  cutoff  other  than  citing  that  it  exists  experimentally.  The 
work  presented  in  this  report  suggests  that  the  drift  of  particles  off  axis 
is  largely  responsible  for  reduced  jet  penetration,  whereby  drifting  particles 
strike  the  target  not  at  the  bottom  of  the  hole  but  somewhere  up  the  hole 
profile.  This  commencement  of  sidewall  impacts  frequently  seems  to  constitute 
a  cutoff  in  jet  penetration,  since  the  overall  depth  of  penetration  is  not 
altered  by  a  sidewall  impact. 

A  method,  therefore,  had  to  be  devised  which  would  allow  penetration  to 
be  a  function  of  drift  velocity.  Three-dimensional  finite  difference  codes 
were  not  considered  due  to  a  coupling  of  availability  and  cost.  Yet,  a  three- 
dimensional  description  of  penetration  was  required  since  the  source  of  altered 
penetration  is  the  spread  of  the  particles  off  of  the  one-dimensional  axis. 
Consequently,  all  aspects  of  particle  penetration  would  have  to  be  considered: 
size,  shape,  and  orientation  of  crater  formation  in  a  three-dimensional  target 
block . 

The  code  developed  in  the  present  work  (referred  to  as  PENJET)  includes 
many  simplifications  and  assumptions.  The  first  assumption  concerns  the 
shape  of  crater  formed  by  an  impinging  jet  particle.  The  chosen  shape  wa 
that  of  an  ellipsoid,  with  the  ellipsoid  being  characterized  by  a  position 
in  three-dimensional  space,  dimensions  along  the  axes,  as  well  as  an  orienta¬ 
tion  in  space  (as  shown  in  Figure  10).  The  position  of  particle  impact  is 
determined  by  using  the  velocity  fan  of  a  round  to  trace  the  particles 
flight  through  three-dimensional  space  and  time  in  order  to  locate  the  particle 
position  with  respect  to  the  instantaneous  hole  profile  (formed  by  earlier 
particles  of  the  jet).  If  a  particle  is  found  to  strike  any  point  within  the 
target,  PENJET  then  proceeds  through  a  series  of  calculations  to  determine 
the  remainder  of  crater  parameters.  Once  the  crater  formed  by  a  particle  is 
effectively  evaluated,  PENJET  proceeds  through  subsequent  particles  of  the 
jet  until  impacts  for  all  digitized  particles  have  been  evaluated.  Thus, 
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P  MAGNITUDE  OF  PARTICLE  PENETRATION 
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Figure  10.  Parameters  Generated  by  PENJET  to  Describe 
Crater  Formed  by  Impinging  Jet  Particle 
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Figure  11.  PENJET  Hole  Profile  Results  from  Ellipsoid  Superposition 
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a  hole  profile  for  a  complete  jet  is  just  the  superposition  of  ellipsoids  from 
each  of  the  particle  impacts.  A  two-dimensional  representation  of  a  three- 
dimensional  hole  profile  is  shown  in  Figure  11  for  a  fictitious  nondrifting 
jet.  When  evaluating  the  crater  formed  by  an  impinging  jet  particle,  crater 
volume  is  a  constraining  factor.  That  is  to  say,  crater  volume  is  calculated 
independently  of  other  crater  parameters  and  is  assumed  to  be  a  function  of 
particle  kinetic  energy  alone.  Subsequent  calculations  for  actual  crater 
dimensions  are,  therefore,  constrained  by  the  allowable  volumetric  target 
displacement  a  given  particle  can  produce.  The  function  governing  the  volume 
energy  relationship  is  a  linear  one,  similar  to  that  adopted  by  DiPer.sio, 
Simon,  and  Merendino. 3  The  distinction  arises  in  the  fact  that  the  DSM  model 
expresses  its  proportionality  constant  as  a  relationship  between  total  jet 
kinetic  energy  and  total  hole  volume.  PENJET, on  the  other  hand,  bases  its 
constant  on  the  energy  of  the  individual  particle.  It  is  not  surprising  to 
find,  therefore,  that  the  DSM  constant  differs  from  its  PENJET  counterpart. 
Consider  also,  that  the  energy  constant  of  PENJET  does  not  strictly  define  an 
exact  amount  of  displaced  target  material.  As  mentioned  previously,  the  hole 
profile  of  PENJET  is  the  superposition  of  ellipsoid  craters  corresponding  to 
each  of  the  jet  particles.  The  orientation  of  these  ellipsoids  will  greatly 
affect  the  spacial  overlap  of  the  ellipsoids,  and  thus  the  overall  volumetric 
target  displacement  of  a  given  jet.  Because  of  this  spacial  overlap  of 
particle  craters,  the  overall  crater  volume  as  calculated  by  PENJET  is  not 
the  summation  of  particle  energies  divided  by  the  energy  constant.  Estimation 
of  the  constant  was  gotten  empirically  through  examination  of  target  blocks 
as  well  as  hole  profile  data. 


The  crater  shape  of  an  individual  impact  is  prescribed  (as  art  ellipsoid) 
and  the  particle's  crater  volume  is  proportional  to  its  own  kinetic  energy. 
Therefore,  knowledge  of  a  particle's  penetration  capability  will  allow  the 
particle's  hole  radius  to  be  computed  and  vice  versa.  In  PENJET,  depth  of 
penetration  is  computed,  with  crater  radius  being  subsequently  determined. 

The  model  of  particle  penetration  used  in  the  code  is  based  on  the  work  of 
Fitzgerald,4  who  draws  on  the  work  of  Kineke5  to  imply  that  crater  formation 
occurs  in  a  direction  perpendicular  to  the  stricken  surface.  Fitzgerald  calls 
this  phenomenon  the  wave  refraction  theory  and  shows  that  it  applied  to  a 
single  particle  striking  a  flat  plate  at  some  obliquity  (Figure  12) .  PENJET 
uses  the  concept  in  a  more  general  sense  by  applying  the  wave  refraction  theory 
to  local  obliquities  as  well.  Local  obliquity  is  defined  as  the  angle  between 
the  particle  flight  axis  and  the  normal  to  the  impact  plane.  Unless  the  im¬ 
pinging  particle  strikes  exactly  at  the  bottom  of  the  hole  profile,  the  re¬ 
fracted  penetration  wave  will  not  be  normal  to  the  original  target  surface 
(Figure  13).  In  the  strict  sense,  the  refracted  wave  theory  should  only  be 
applicable  in  the  regime  of  hydrodynamic  penetration.  A  large  majority  of 
narticle  impacts  studied  did  not  fit  the  qualification  for  hydrodynamic  pene¬ 
tration.  For  these  cases,  the  direction  of  penetration  would  be  somewhere 
between  the  particle  flight  axis  and  the  refracted  wave  axis.  It  was  later 
noted  that  the  penetration  of  a  particle  typically  varied  less  than  5%  because 
of  the  refracted  wave  assumption. 

4£.  Fitzgerald ,  Particle  Waves  cmd  Deformation  in  Crystalline  Solids .  Inter- 
science  Publishers,  New  York,  1966, 

V.  Kineke,  Jr.  "An  Experimental  Study  of  Crater  Formation  in  Metallic  Tar- 
gets  ."Proa.  4th  Symposium  Hypervelocitu  Impact.  Vol.  I  (Held  Eqlin  AFB, 
Florvda,  April  26-28,  1960). 
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Figure  12.  Crater  Formation  According  to  Refracted  Particle  Wave  Theory 


EFFECTIVE  PARTICLE 
LENGTH  (Leff) 


Figure  14.  Concept  of  Effective  Length  as  Seen  by  Impact  Plane 


Jet  penetration  in  the  continuous  mode  is  quite  a  different  matter  from 
in  the  particulate  mode.  Modifications,  therefore,  had  to  be  included  to 
describe  continuous  jet  penetration.  The  criterion  on  which  to  base  mode  of 
penetration  is  that  of  breakup  time.  PENJET  uses  a  linear  breakup  equation 
whereby  the  jet  tip  particulates  at  some  t^n  and  subsequent  particle  breakup 

occurs  at  regular  At  intervals.  Estimates  of  jet  particulation  are  obtained 

through  examination  of  the  radiographs  for  the  three  different  exposures. 

These  estimates  are  then  used  in  a  least  squares  fit  to  evaluate  t  .  and  At. 

^  mm 

Based  on  the  estimated  breakup  time,  it  can  be  established  whether  a  given 
particle  will  be  penetrating  continuously  or  particulately  at  a  given  standoff. 
If  continuous,  the  length  of  the  particle  is  re-evaluated  such  that: 
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IMPACT 
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Also,  it  is  assumed  that  jet  penetrating  in  the  continuous  mode  does  not 
experience  the  refracted  particle  wave  of  Figure  12.  This  must,  of  course  be 
the  case  since  the  continuous  particle  is  essentially  beginning  its  penetra¬ 
tion  exactly  when  and  where  the  previous  particle  leaves  off. 

PF.NJET  also  considers  particle  orientation  with  respect  to  the  stricken 
crater  surface.  The  assumed  shape  of  a  shaped-charge  particle  is  cylindrical, 
being  characterized  by  a  length  and  diameter.  As  local  obliquity  is  increased, 
the  length  of  the  particle  as  seen  by  the  impact  point  is  effectively  changed. 
(See  Figure  14.)  The  extreme  case  is  when  local  obliquity  approaches  90°. 

In  this  case,  the  effective  length  of  the  particle  approaches  the  diameter  of 
the  particle.  Hydrodynamically,  the  penetration  of  a  particle  is  proportional 
to  the  particle  length.  Thus,  the  effective  particle  length  replaces  the 
particle  length  in  penetration  calculations. 

The  presence  of  obliquity  raises  more  questions.  The  larger  the  obliquity, 
the  smaller  the  component  of  particle  velocity  perpendicular  to  the  stricken 
surface.  Once  in  the  hypervelocity  regime,  large  changes  in  particle  velocity 
produce  relatively  small  changes  in  penetration.  If  particle  velocity  is  too 
small,  however,  the  strength  of  the  target  and  jet  play  a  major  role  in  reduc¬ 
ing  the  penetration  of  that  particle.  Besides  the  effective  particle  length 
that  the  target  "sees",  it  is  believed  that  an  effective  velocity  also  plays 
a  role  in  determining  depth  of  penetration.  Whereas  penetration  is  assumed  to 
be  directly  proportional  to  effective  length,  the  role  that  effective  velocity 
plays  is  not  nearly  so  profound.  As  long  as  the  effective  velocity  is  in 
the  hypervelocity  regime  (arbitrarily  chosen  as  any  velocity  above  2mm/psec), 
penetration  is  assumed  greater  than  the  effective  particle  length  for  copper 
jet  into  an  RHA  target.  As  effective  velocity  increases,  penetration  is 
assumed  to  asymptotically  approach  the  hydrodynamical ly  idealized  value. 

That  is:  ,  ,  ,  , - ; - 

AP  -+  where  y  =  /  P^/  P •. 


PENETRATION/PARTICLE 


1.5  i- 


O 

Z 


1.0 


0.5  h 


0.0 


r 


P/L=1.0- 


(Veff-2.0) 


2  ' 


s'*'"  HYDRODYNAMIC  + 


STRENGTH  TERMS 


0.72 


I 


I  (Vcff  "0.8) 2 
/  P/L  =  - 


P/L  =  0.0  / 

I— L  l.J-1, 


0.72 
I -L  I 


1—1  J.  — L  „  1—  I  i  LLLIJ 


12  3  4 

EFFECTIVE  VELOCITY  (mm/^s) 


Figure  15.  The  Role  of  Effective  Velocity  in  Penetration  Calculations 
Involving  Copper  Jets  Against  Steel  Targets 
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The  asymptotic  approach  is  generated  using  the  hydrodynamic  equations,  includ¬ 
ing  strength  terms.  For  effective  velocities  below  2mm/usec,  penetration 
rapidly  decreases  in  a  manner  expressed  in  Figure  15  so  that  no  penetration 
is  induced  when  the  effective  particle  velocity  is  lowered  beyond  .8  mm/ usee. 

Data  for  the  effective  velocity  curve  in  the  low  velocity  regime  was  gathered 
from  work  by  Weirauch.6  Impacts  of  this  nature  may  be  the  source  of  particle 
residue  shown  by  Simon  and  PiPersio7  to  accumulate  at  the  bottom  of  the 
hole  profile.  PF.N JET, however,  does  not  consider  this  type  of  interference. 

VI.  APPLYING  THE  CODE  TO  DATA 

Given  a  set  of  digitized  radiographic  data,  PENJET  will  go  through  all 
of  the  mentioned  calculations  for  each  digitized  particle  of  the  jet  to  even¬ 
tually  create  a  three-dimensional  hole  profile  for  a  given  round  at  a  given 
standoff.  If  the  code  is  told  to  use  the  data  of  a  given  round  to  perform 
penetration  calculations  at  a  variety  of  standoffs,  a  theoretical  penetration 
standoff  curve  can  be  generated  for  each  round.  This  curve  is,  of  course,  only 
valid  for  a  particular  round.  Thus,  each  generated  curve  is  as  individual  as 
the  data  from  which  it  was  generated,  and  does  not,  in  general,  pertain  to  all 
rounds  of  the  same  type. 

The  code  has  both  its  good  and  bad  aspects.  A  major  consequence  is  that 
tilt  has  no  direct  effect  on  the  code's  ability  to  calculate  penetration.  This 
happens  because  tilt  does  not  misalign  the  particles.  Rather,  it  only  redefines  the 
charge  axis.  The  particles  still  drift  with  respect  to  each  other  in  the  same 
fashion  despite  measurements  including  tilt.  Much  emphasis  was  previously 
placed  on  tilt  for  the  reason  that  it  would  be  advantageous  to  obtain  drift 
distributions  unbiased  by  tilt.  In  terms  of  the  penetration  calculations, 
though,  tilt  of  and  by  itself  has  no  effect. 

Tilt  does,  however,  indirectly  affect  the  situation  in  terms  of  reference 
drift  calculations.  As  relative  drift  velocity  (as  calculated  from  the  radio¬ 
graphs)  is  the  superposition  of  tilted  and  actual  drift,  abnormal  relative 
drift  distributions  can  occasionally  result  depending  on  how  the  two  components 
of  measured  drift  interact  with  each  other.  Correspondingly,  abnormal  relative 
drift  distributions  lead  to  abnormal  extrapolations  of  reference  drift,  thus 
making  the  calculated  reference  drift  suspect.  After  all,  the  calculation  of 
reference  drift  is  only  an  educated  approximation:  Direction  of  the  drift  is  based 
an  jet  particle  orientation,  while  magnitude  of  drift  velocity  is  based  on  Equa¬ 
tion  (9).  Therefore,  unlike  relative  drift  computations,  which  are  as  good  as 
accuracy  of  measurement,  reference  drift  is  only  as  good  as  the  governing 
assumptions.  Thus,  the  poorer  the  approximation  at  reference  drift,  tbe  poorer  the 
quality  of  penetration  calculations. 

6  G.  eihvauch ,  "The.  Behavior  of  Copper  Pins  Upon  Impacting  Various  Materials 
With  Velocities  Between  50  m/s  and  1650  m/s>"  (Doctoral  Thesis)  University 
of  Karlsruhe ,  1971. 

7  J.  Simon,  R.  DiPersio,  " Experimental  Verification  of  Standoff  Effects  on 
Shaped-Charge  Jet  Cutoff  in  Solid  Targets Ballistic  Research  Laboratory 
Memorandum  No.  1976,  May  1969  (AD  854396). 


Another  drawback  to  the  code  results  from  the  radiographs  themselves. 

The  code  can  only  analyze  particles  that  have  been  digitized.  In  ">r»s*  radio¬ 
graphs  some  portion  of  the  jet  that  is  visible  has  not  yet  particulated. 
Attempting  to  separate  the  continuous  jet  into  digitizable  portions  has  met 
with  some  success,  though  likely  at  the  cost  of  accuracy  of  particle  length, 
since  the  continuous  particles  are  still  stretching  when  digitized.  Also, 
most  radiographs  don't  even  contain  the  complete  jet,  for  there  is  still  jet 
matter  that  hasn't  reached  the  plane  of  the  film  by  the  time  it  is  exposed. 
Thus,  there  is  no  way  to  digitize  these  particles.  For  the  rounds  studied,  it 
is  estimated  that  over  40%  of  the  jet  was  routinely  nondigitizable. 

To  what  extent  the  ability  to  digitize  particles  affects  the  penetration 
code  predictions  depends  on  the  quality  of  the  jet.  In  fact,  the  poorer  the 
jet  is,  the  less  predictions  are  affected.  If  the  undigitized  particles  drift 
far  enough  off  course  from  the  front  portion  of  the  jet,  they  will  never  make 
it  to  the  bottom  of  the  hole  profile  to  add  to  penetration.  In  this  situa¬ 
tion,  calculations  based  only  on  the  particles  digitized  should  produce 
results  compatible  with  calculations  involving  all  particles.  Since  distance 
of  deviation  is  directly  proportional  to  time  in  flight  which  is  directly 
proportional  to  standoff,  longer  standoffs  will  also  decrease  the  penetration 
effectiveness  of  the  rear  of  the  jet.  Thus,  penetration  predictions  at  short 
standoff  are  subject  to  unavoidable  error  since  the  undigitized  jet  portion 
is  still  capable  of  increasing  jet  penetration.  The  fact  that  the  penetration 
effectiveness  of  the  jet  rear  is  predicted  to  decrease  with  increasing 
standoff  correlates  with  the  DSM  concept  of  minimum  jet  velocity  increasing 
with  standoff.  3 

Having  kept  these  reservations  in  mind,  the  results  of  the  code  based  on 
data  from  four  rounds  are  illustrated  in  Figures  16-19.  The  only  factor  govern¬ 
ing  the  choice  of  these  particular  four  rounds  was  penetration;  that  is, 
rounds  were  chosen  with  varying  penetrations  ranging  from  poor  to  good.  The 
code  was  run  at  standoff  increments  out  to  38  CD.  The  actual  rounds  were 
fired  at  24  CD  standoff  to  allow  radiography,  and  the  actual  penetrations 
achieved  are  shown  on  each  of  the  graphs  as  a  single  datum.  Steplike  varia¬ 
tions  in  the  code's  predictions  result  because  of  the  discreteness  of  the 
impacting  particles:  Either  a  particle  made  it  to  the  bottom  of  a  hole  profile 
or  it  did  not.  Major  fluctuations  in  the  code  predictions  can  be  taken  to 
imply  uncertainty  in  the  jet  performance.  Uncertainty  as  to  whether  two  or 
three  given  particles  make  it  to  the  bottom  of  the  hole  profile  can  create 
penetration  variations  in  excess  of  1/2  CD.  A  curve  was  drawn  to  fit  the 
code's  predictions  on  each  of  the  rounds,  though  a  band  width  of  uncertainty 
should  be  realized  as  well.  This  curve  is  PENJET's  prediction  of  penetration, 
assuming  that  one  could  fire  the  exact  same  round  over  and  over  again  at 
various  standoffs.  Since  we  know  that  the  predictions  at  short  standoff  are 
inaccurate  due  to  nondigitizable  jet,  the  code's  predictions  on  Figures  lb- 
19  are  generally  unreliable  below  10  CD.  The  curve,  however,  should  be  valid 
at  the  standoff  of  24  CD  at  which  the  rounds  wete  actually  fired. 
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Figure  16.  Analysis  of  Round  2331 
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Figure  17.  Analysis  of  Round  2332 


Figure  18.  Analysis  of  Round  2334 
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Figure  19.  Analysis  of  Round  2336 


Figure  20.  Routine  Performance  for  Type  of  Round  Analyzed 


All  of  the  codes 's  predictions  were  in  reasonable  agreement  with  experi¬ 
ment  except  for  Round  2332.  Upon  studying  this  round  for  possible  causes  of 
discrepancy,  it  was  noted  that  the  round  uniquely  had  what  shall  be  called 
a  switchback  angular  spiral.  This  falls  into  the  category  of  abnormal  rela¬ 
tive  drift  distributions.  Most  jets  contain  drift  distributions  such  that 
the  jet  particles  spiral  in  only  one  direction  (either  clockwise  or  counter¬ 
clockwise)  throughout  the  jet.  Figure  5,  for  example,  shows  a  jet  spiralling 
counterclockwise  from  jet  tip  to  tail.  A  switchback  angular  spiral  is  one 
in  which  one  part  of  the  jet  spirals  in  one  direction,  while  the  subsequent 
portion  of  the  jet  spirals  in  the  opposite  direction.  Switchback  is  thought 
to  be  a  product  of  abnormal  tilt/drift  combinations  such  that  subsequent 
computation  of  reference  drift  might  be  quite  inaccurate.  As  set  forth  before, 
an  inaccurate  reference  extrapolation  leads  to  inaccurate  penetration  calcula¬ 
tions. 

Figure  20  shows  how  the  actual  weapon  routinely  performs.  In  order  to 
test  whether  PENJET  could  predict  the  average  penetration  at  standoffs  near 
the  peak  performance  of  the  round,  a  radiograph  containing  a  complete  and 
fully  particulated  jet  would  be  required  for  digitization.  As  mentioned  pre¬ 
viously,  no  available  radiographs  fulfilled  these  requirements.  A  radiograph 
was  obtained,  however,  of  a  complete  jet  which  was  not  fully  particulated. 
Though  not  completely  digitizable,  the  presence  of  a  complete  jet  greatly 
facilitated  estimation  of  many  jet  parameters  (e.g.  velocity  of  jet  rear, 
drift  with  respect  to  jet  rear,  and  total  jet  length).  By  this  means,  the 
rear  of  the  jet  which  was  not  digitizable  could  be  simulated  into  particles 
based  on  analysis  of  the  radiograph  containing  the  continuous  jet  rear. 

Figure  21  shows  the  PENJET  prediction  of  the  round  that  includes  the  jet  rear 
simulation.  The  actual  datum  as  well  as  optimum  performance  for  this  type  of 
round  are  included  in  order  to  put  PENJET  calculations  in  perspective. 

The  importance  of  a  complete  jet  at  short  standoff  is  shown  in  Figure  22. 
The  penetration  standoff  curve  for  Round  2331  is  shown  here.  Superimposed 
upon  it  are  curves  showing  the  amount  of  residual  penetration  resulting  from 
different  fractions  of  the  complete  jet.  The  points  at  which  the  curves  meet 
signify  that  the  number  of  particles  producing  the  superimposed  curves  are 
also  the  same  number  of  particles  producing  the  total  residual  penetration  at 
that  point.  For  example,  PENJET' s  predictions  at  11.5  CD  standoff  imply  that 
predicting  penetration  using  40  particles  would  be  just  as  accurate  as  using 
51  particles.  Anything  less  than  40  particles  at  that  standoff  would,  however, 
produce  inaccurate  results.  For  the  24  CD  standoff  at  which  the  round  was 
actually  fired,  less  than  20  digitized  particles  would  be  required  for  an 
accurate  description  of  penetration.  At  short  standoff,  however,  even  51 
particles  do  not  completely  describe  the  penetration  capabilities  of  Round 
2331.  Estimates  are  that  below  8  CD  standoff,  particles  beyond  particle  51 
were  also  capable  of  adding  to  penetration.  For  better  quality  rounds,  though, 
as  illustrated  by  Round  2937  in  Figure  23,  more  particles  are  required  even 
at  longer  standoffs  to  accurately  describe  penetration.  For  example,  more 
than  50  digitized  particles  would  be  required  for  standoffs  out  to  25  CD  for 
this  particular  round. 
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Figure  21.  Analysis  of  Round  2937  with  the  Included  Jet  Rear  Simulation 
Note:  Optimum  Performance  Data  Transferred  from  Figure  20. 


Figure  23.  Penetration  Capability  of  Distinct  Jet  Portions  of  Round  293? 


Figure  25. 


PENJET  Predictions  of  Hole  Profile  for  Round  2937  at 
CD  Standoff.  Penetration  in  Cone  Diameters. 


23 


43 


IMPACT  DEPTH  (CD) 


PARTICLE  NUMBER 


Figure  26. 
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Figure  27.  Orthogonal  Views  of  Projected  Hole  Profile  of  Round  2331  at  18  CD 
Standoff  as  Predicted  by  PENJET.  Penetration  in  Cone  Diameters. 
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Figure  31.  Orthogonal  Views  of  Projected  Hole  Profile  of  Round  2331  at  38  CD 
Standoff  as  Predicted  by  PENJET.  Penetration  in  Cone  Diameters. 


Besides  showing  a  semblance  of  agreement  between  code  and  experiment. 
Figures  16-19  show  that  drift  alone  can  easily  account  for  rapid  decreases  in 
penetration  as  is  often  experienced  by  rounds  at  extended  standoffs.  The 
process  by  which  penetration  is  decreased  is  one  of  sidewall  impacts.  This 
process  was  briefly  mentioned  before  as  a  phenomenon  which  sometimes  renders 
a  particle  ineffectual.  The  majority  of  sidewall  impacts,  though,  are  not  so 
oblique  as  to  disrupt  penetration;  most  sidewall  impacts  go  into  creating 
wider  hole  profiles. 

As  an  example  of  the  devastating  effects  of  drift  velocity  on  jet  penetra¬ 
tion,  a  comparison  of  PENJET's  hole  profile  predictions  is  shown  for  two 
separate  rounds  in  Figures  24  and  25.  The  plots  represent  orthogonal  views 
of  the  target  block  as  described  by  PENJET.  Both  sets  of  calculations  were 
performed  at  the  same  standoff  (23  CD).  The  stark  contrast  of  jet  perform¬ 
ance  is  noted  not  only  in  terms  of  depth  of  penetration,  but  also  in  terms  of 
crater  width.  Sidewall  impacts  do  not  occur  until  much  later  in  the  penetration 
process  for  Round  2937,  thus  allowing  more  of  the  jet  to  contribute  to  the 
overall  penetration. 

Figure  26  is  PENJET's  prediction  as  to  where  in  the  hole  profile  each 
particle  of  Round  2331  strikes  the  target  block  at  18  CD  standoff.  Note  that 
particle  number  is  in  order  of  formation.  (Jet  tip  is  particle  1.)  It  is  seen 
that  once  a  particle  is  blocked  from  reaching  the  bottom  of  the  hole,  it  will 
usually  begin  the  widening  process  somewhere  up  the  profile.  Subsequent 
particles  try  to  work  their  way  back  down  again.  This  process  continues  until 
no  jet  is  left  with  which  to  penetrate  the  target.  A  plot  of  the  hole  profile 
as  generated  by  PENJET  is  shown  in  Figure  27  as  two  orthogonal  views  of  the 
target  block.  Graphs  and  plots  for  the  same  round  are  shown  again  at  longer 
standoffs  in  Figures  28-31.  It  is  seen  that  penetration  is  progressively 
lower  because  sidewall  impacts  occur  earlier  in  the  penetration  process. 
Consequently,  crater  diameter  progressively  increases  at  longer  standoffs. 

Thus,  for  Round  2331  at  a  standoff  of  23  CD  the  total  penetration  capability 
of  the  jet  as  predicted  by  PENJET  is  contained  within  the  first  18  particles. 

The  profile  at  38  CD  sharply  illustrates  the  effects  of  drift  velocity 
on  jet  penetration.  The  spread  of  the  particles  was  so  great  as  to  cause  the 
beginnings  of  a  dual  crater  within  1  CD  of  the  target  surface.  Thus,  the 
ineffectiveness  of  the  jet  rear  is  seen.  Any  nondigitized  portion  of 
Round  2331  would  probably  produce  its  penetration  in  the  vicinity  of  the 
particle  51  crater,  far  from  the  point  of  maximum  penetration. 

Thus,  the  decreases  in  penetration  as  predicted  by  PENJET  are  due  mainly 
to  the  presence  of  drift  velocities  among  the  particles.  Drift  causes  scatter¬ 
ing  of  the  particles  and  eventually,  at  some  critical  standoff,  elimination 
of  the  particles  from  the  crater  deepening  process  via  the  mechanism  of  side- 
wall  impacts. 
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VII.  CONCLUSIONS 


The  ability  to  compute  drift  velocities  for  the  particles  of  a  shaped  - 
charge  jet  is  thought  to  be  significant.  In  its  own  right,  quantification  of 
drift  velocity  aids  in  understanding  the  nature  of  the  nonideal  jet.  A 
degree  of  imperfection  may  also  be  placed  upon  the  jet  in  hopes  of  classifying 
a  jet  with  regard  to  its  quality. 

The  concept  of  how  a  jet  penetrates  is  better  understood  when  drift 
velocity  is  considered  in  the  penetration  picture.  The  process  of  sidewall 
impacts  as  a  hinderancc  to  penetration  can  be  verified  by  inspection  of  actual 
target  blocks.  This  process  is  reminiscent  of  the  jet  cutoff  velocity  which 
varies  with  standoff  as  proposed  by  Dipersio,  Simon,  and  Merendino,3  in  that 
the  velocity  of  the  particle  striking  deepest  in  the  hole  profile  can  be  con¬ 
sidered  a  cutoff  velocity  with  respect  to  jet  penetration.  Implications  here 
are  that  a  jet  could  approach  ideal  penetration  had  drift  velocity  not  been 
present.  Larger  craters  would  also  make  jets  more  insensitive  to  drift. 

Weaker  targets  allow  for  a  wider  crater  for  a  given  jet  particle.  Thus,  the 
reduced  probability  of  sidewall  impact  would  allow  greater  penetration. 

Compared  to  armor,  a  target  of  mild  steel  hydrodynamically  should  not 
affect  the  penetration  of  a  given  jet.  It  will  affect  penetration, however , 
when  drift  velocity  is  present  in  the  jet.  DiPersio  and  Simon  8  verified  this 
experimentally  by  firing  jets  into  steel  targets  of  varying  hardness.  Increas¬ 
ing  crater  size  could  alternately  be  accomplished  by  means  of  increasing  the 
jet's  energy.  If  dynamic  impact  did  not  induce  additional  drift  in  the  jet, 
projecting  the  round  into  the  target;  as  is  the  case  under  field  conditions  , 
should  produce  greater  penetration  than  the  static  detonation  of  the  same 
round.  This  cannot  be  verified  experimentally,  however,  due  mainly  to  the 
complexity  of  dynamic  impact.  Besides  possessing  possible  yaw  and  pitch,  the 
dynamic  round  is  moving  with  respect  to  the  target  during  the  jet  formation 
nroccss.  This  relative  motion  can  manifest  itself  in  the  form  of  drift. 
Dynamically  detonated  rounds  do,  however,  create  wider  craters  than  their 
static  counterparts  as  expected. 

It  is  suggested  that  future  investigations  attempt  to  analyze  drift  from 
the  viewpoint  of  warhead  causes.  The  first  step  in  understanding  penetration 
in  terms  of  drift  has  hopefully  been  set  forth  in  this  report.  The  ability  to 
analyze  drift  as  a  function  of  the  warhead  would  be  the  link  between  the 
imperfect  warhead  and  nonideal  penetration.  Eventually,  the  possibility 
exists  where  warhead  tolerances  may  be  used  to  place  limits  upon  the  drift 
velocity  distributions  for  a  round  so  that  penetration  might  be  expressible 
as  a  function  of  charge  tolerances. 


6R.  DiPersio ,  J.  Simon ,  "The  Effect  of  Target  Hardness  on  the  Penetration  Capa¬ 
bility  of  Shaped -Charge  Jets,"  BRL  Report  No.  1408 ,  July  1968  (AD  838991). 
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The  drift  velocity  and  PENJET  programs  were  written  as  direct  additions 
to  Blische's  and  Simons'  radiograph-data-reduction  program  (Reference  1) 
introduced  in  BRL  Report  2330.  Therefore,  a  listing  of  PENJET  must 
be  accompanied  by  the  program  which  generates  data  for  use  within  PENJET.  The 
right-hand  edge  of  the  listing  reveals  where  the  statements  in  question 
originated.  "GETREADY"  are  nonexecutable  statements  required  for  BRL  Report 
2330,  PENJET,  or  both.  "BRL  2330"  is  the  bulk  of  Blische's  and  Simmons' 
program  with  minor  revisions  for  PENJET  applicability.  "JETBREAK"  is  a 
revision  for  jet  breakup  calculations.  "DRIFT"  is  used  to  calculate  the  drift 
distribution  for  a  round  while  "PENJET"  is  the  jet  penetration  program.  The 
program  is  written  in  FORTRAN  IV.  An  explanation  of  the  program  variables  is 
inserted  after  the  program  listing. 

It  should  be  noted  that  even  though  Blische's  and  Simmons'  program  is 
versatile  enough  to  run  using  data  from  only  two  X-ray  flash  exposures, 

PENJET  is  geared  to  the  standard  3  flash  setup  and  will  require  modifications 
to  run  otherwise. 

Explanation  for  the  program  variable  names  found  in  the  program  segment 
identified  as  "BRL  2330"  can  be  found  in  BRL  Report  2330  (Reference  1) . 

Should  one  wish  to  identify  variables  within  the  program  segment  "GETREADY," 
the  other  segments  should  be  consulted,  as  all  "GETREADY"  variables  are  con¬ 
tained  within  other  segments.  All  other  significant  variables  are  listed 
under  their  appropriate  program  segment. 
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TROGRAH  MAIN  (  INPUT, OUTPUT*  TAPES  .INPUT,  TAP  16*  OUTPUT  ) 
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OIHENSION  AA(2,3>,CCt2>*AR 13001, AT  13001. S1«I 2) ,TT tZ 3 
OIHENSION  OS  (3).TH£TAI100).yETRAV 1 100). Til 1001 *12(100). Tit  1001* 
t  REFI3) 

OIHENSION  AVSTtl00)*BETAI100) 

Of PENSION  PHI(100)*PHIR(]00)*VTABSI100) 

01 HENS  ION  EPS  I  6*100) *  VTI 6) 

01HENS10N  VX(100),VYI100> 

07NENS10N  HI?) 

OIHENSION  VT8I100) 

01  HENS  ION  7llOO)*X(100),Y(100)*Z(100),U(100)*RADtUS(100)*P€N(100)* 
1  HV0LI100) 

OIHENSION  NCCdOO) 

01  HENS  ION  ZSUB(100),PENT(100t,PENNAXtl00) 

OIHENSION  Vt  100*3>*VimOO),VAI100,ll*VARtlOO>tVPI100,3l,VPN(100t  • 

1  RGRAOI 100*  3 )*RGH(  100  )«GAAD(100,3I*AR(100*3)*  PR  (100*  3) 

OIHENSION  Rl(100)*RVELI100> 

OIHENSION  OU100)*KWAU100)«TAU(100) 

OIHENSION  TP  ( 100  )-*RD  1 100) 

REAL  KOI?* KD23.KD13 
PE  At  I  NCR 

OIHENSION  SA(10),SBI10).SC(10).SO(10>,SE(10>*$FI10>.S11(10I» 

1SL2 1 10). St  31 10). ST (10), ll( 100), 12(100). 13(100) 

OIHENSION  PI100)*T0TP(  100)  *S(JN01A  ( 100)  ,  SUNIENI  100)*XDf  A(  100)>0EtV( 
1 100)*SDELV(100)*  SUNOEt (100) 

D I  HE  NS  1  ON  XV0K3.100) 

OIHENSION  AK100.100) 

DATA  SA(1)*SA(2)*SA(3|/1 OH VELOCITY  I • 10HHN/H ICROSE, 3HC )»/ 

DATA  SB(1)*SB(2)*S8( 3) I10HCUNUL AT IVE * 10H  NASS  <GRA,AMNS)»/ 

DATA  SC(1).SC(2).SC(3)*SC(A) /10HP0S1T ION  A.10HL0NG  JET  1.1CHENGTH 
1(HN).1H>/ 

DATA  SO(l).SO(2).S0(3)/10HCUNULATIVE.10H  K.E.  ( JOU*  3HLES I > / 

DATA  SE(1)*SE(2).SE(3)*SE(A) 710HD1ST ANCE  E.10HR0H  CHARGE. 1 CH  BASE 
1 (HH),1H>/ 

DATA  3F(1),SF(2)*SF(3) /l OH BREAK -UP  T.IOHIRC  (HICRO, 3HSEC )» / 

DATA  SL1 ( 1 ) / BHEL ASH  1»/ 

DATA  St  2 ( 1  )  / BHF  L  ASH  2 »/ 

DATA  St 3 ( 1 ) /BHF L  ASH  3»/ 

DATA  HPA,HPAA,HPART,HCO.NCON»HCONT/AHPART,AHtCUt,AHATED.AH  CON.AHT 
1 INU, AHOUS  / 

1  FORHAT(3I3*FS.2*6F10.S) 

2  FORNAT(YFIO.S) 

6  FORNAT(BFIO.S) 

21  FORHAT (215) 

30  F0BHAT(1?13,ZX,I1,2X,313) 

RE  AD( B*  2 1 )  ICASES*  tCL 
IF  (  EOF  ( 3)  )  23.23 
23  00  BOO  IJ-l.TCASES 

2?  RE AO (3. I INROUND. NPART .NFLASH.RHC* (XHAG( I ) ,!• I , 3) .Ft ASHl.Fl ASH2, 
IFLASH3 

IF  I  EOF (9 ) )  2A.2A 
2A  ENCQ0E(21.2O*ST(D)  NROUNO 
20  FORNATaOHROUNO  NUNB.AMER  ,  IS.2H  >1 

RE  ADI  5. 2)  F0CUS1. FOCUS?* FOCUS 3»F1A*F2A*F3A.FA A 
IF  (EOF  (3) )  2  3, 2 B 

23  READ) 3,6)  F1B,F2B,F3B.FA8,F1C*F?C,F3C,FAC 
IF (EOF (3 ) )  26,26 
26  RO.RHO 

<IFIICI.NE.1I  GO  TO  29 


CETREADY 

2 

GETREAOY 

3 

getreaoy 

A 

GETREAOY 

3 

GETREAOY 

6 

getreaoy 

r 

GETREAOY 

6 

GETREAOY 

9 

GETREAOY 

10 

GETREAOY 

11 

GETREAOY 

12 

GETREAOY 

13 

GETREAOY 

1A 

GETREAOY 

13 

GETREAOY 

16 

GETREAOY 

IT 

GETREAOY 

IB 

GETREAOY 

19 

GETREAOY 

20 

GETREAOY 

21 

GETREAOY 

22 

GETRE AOY 

23 

GETREAOY 

2A 

GETREAOY 

23 

GETREAOY 

26 

GETREAOY 

27 

GETREAOY 

26 

GETREAOY 

29 

GETREAOY 

30 

G'TREADY 

31 

GETREAOY 

32 

GETREADY 

33 

getreaoy 

3A 

GETREAOY 

33 

GETREAOY 

36 

GETREAOY 

37 

getreaoy 

36 

GETREAOY 

39 

GETREADY 

AO 

GETREAOY 

A1 

GETREADY 

A2 

GFT  READY 

A3 

GETREADY 

AA 

GETREAOY 

A3 

BRL2330 

2 

BR12330 

3 

BR12330 

A 

BR12330 

5 

BR12330 

6 

BR17330 

7 

BRl 2330 

B 

BRL  2330 

9 

BRL  2330 

10 

BR12330 

11 

BRL 2330 

12 

BR12330 

13 

BRl 2330 

1A 

BRl 2330 

13 

BR12330 

16 

BRL  2330 

17 

BRL  2330 

IB 

BRL2330 

19 

BRl 2330 

20 

58 


ROUND  INFORMATION  CARO  NEXT 

29  PRINT  5 

PRINT  3,  NROUND.  XHO, ( XHAG 1 N l,N. 1, 3 ), FOCUS  1. FOCUS2. FOCUS  3, F l ASH1 , FL 
1ASH2.FLASH3 

3  FORMAT! !/#//• 20X, 'ROUND  NUMBER  •»U»  //»  20X, 'LINER  OENSITYIGM/CC  »- 
1  ’  »FA,  1./.20X,  'MAGNIFICATION  F  AC  TOR-'  •  SF  « .3.  /  ,  20X, ’DISTANCE  FROM  L 
2TNER  BASE  TO  FOCAL  POINTIMM! • > /< 2SX, • Ft  ASH  1-  • , Ft. I. S 2fX, 'FI  ASH 
32-  SF6.1,/,29X, 'FLASH  3-  '. F6. I. /.20X ,' DELAY  TIMES  (MICROSEC )• t /• 
A25X, ’FLASH  1-  '.F6.1,S25X-,’FIASH  2-  SF6 .1,/, 2JX, 'FLASH  3-  SF6.1 
31 

RHO’RO 
PRINT  3 
3  F  DR  MAT ( 1  HI ) 

PRINT  A 

A  FOR HATI21X, 'PARTICLE  AVG.  VELOCITY  TOTAL  JET  BREAK-UP' »J»2 

12X, 'NUMBER SAX,’  (MN/H1CROSEC I S 3X, 'LENGTH! HN1 S 3X.’ I PICROSEC > • , / I 
PI-3.1A139 
00  130  I-1,NFLASH 
DO  1A9  J>1,NPART 

READ(9,30I!IXIK)-,IV!KSK«l*6),l,IFTLH,  I  ROUND#  IPART 
IF  I  EOF ( 31 >  32,32 

32  DO  33  M«li6 

Z XI H>-F10AT(TXIH Ilf  13.9906 

33  ZVIMI«FLOAT(IVIHI 1/13*3906 
Alll,|)>!  IZYt2l»ZY(3ll/2.l*XMACill 
1FIZXII1I.E0. 2XI2II  GO  TO  AO 

Rl>.3*SaRT((2X(II-2X(6ll**2«(rr<n-2r(6t>*P2l*XNA6m 
R2-.5»SORT((IX(2t-IX!3))*P2,|2T(2)-2Y(5) ) **2 ) *XHAG(  I  I 
R3-.9*SORTI!ZX(3l-ZX(AI1**2«(ZY(3l-ZVIAII*P2IPXNA6(II 
P1X*IZXI1)*ZXI6) 1/2. 

P1Y-(ZY(1I«ZY(61 1/2. 

P2X.(ZX(3)'ZX(AI 1/2. 

P2VMZY(3I*ZV(AI  1/2. 

P3XMZX(2l’ZX(9ll/2. 

P3Y>(ZYI2t*ZYIf  11/2. 

XHl«SQRT(IPlX-P3XI**2»(PlV-P3YI**2t*XRACttt 
XM2-S0RT! !P3X-P2X)P*2a(P3Y-P2Y)PP2)PXMAGII> 

P1Z-P1X 
P1X-P2X 
60T0  A3 

ao  Pix-zxm 

A3  1F!L.EQ.2I  GOTO  33 
IFIL.E0.3I  GOTO  63 
IFdFILM.E0.2l  60  TO  AT 
IP! IFTLM.EQ.3 I  GO  TO  AB 
IF  I  IF 1LN.EQ. A)  GO  TO  A9 
SI! J)-F0CUS1-IFIA-PIX)PXMAG(II 
GO  TO  13 

AT  III JI«F0CUS1-(F2A-P1XI*XHA6(II 
GO  TO  If 

AS  SlIJI’FQCUSl’IF 3A,P1XI*XHAG(  1 1 
GO  TO  13 

A9  S1(IM«F0CUS1*(FAA»P1X>*XHAG(  I) 

13  IF!2mi.C0.2X(2M  60  TO  30 
LIIJI’O 
GOTO  T3 

30  L1(JI*1 
GOTO  BO 

33  1FIIFUN.E0.2)  60  TO  fT 
ITIIFUM.E0.3I  GO  TO  SB 
IF (TFILM.EQ. A  I  60  TO  39 
$2(JI’F0CUS2-!F1B-P1XI*XHAGI1I 


BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BR  L  23  30 
BRL2330 
BRL2330 
BRL2330 
6RL2330 
BRL2330 
9X12330 
BRL2330 
9X12330 
BRL2330 
9X12330 
BRL2330 
BR  L  23  30 
BRL2330 
9X12330 
BRL2330 
9X12330 
BRL2330 
9X12330 
BRL2330 
9X12330 
BRL2330 
9RL2330 
BRL2330 
BRL2330 
9X12330 
BRL2330 
BRL2330 
BRL2330 
9RL2330 
BRL2330 
BRL2330 
6RL2330 
9X12330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
9X12330 
9X12330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
9X1 2330 
BR12330 
9RL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
BRL2330 
9X12330 
BRL2330 


CO  TO  1C 

*7  S2» JUFOCUS2-1F II-MDMIUdl) 

60  TO  1C 

60  52<J>>F0CUS2t(Fl(«»lX)*XHA6m 
GO  TO  1C 

S«  S2m«FOCUS2*(F*6*FlX)*XNACII) 

1C  VII  J)>(S2<.l)-Slt6))MFLASH2-FlASHl> 

1FI2XI1). 20.20(21)  CO  TO  CO 
121  J».0 
GOTO  75 
60  12IJN1 
GOTO  60 

66  1TINFIASN.I0.2I  GOTO  60 
IF(|FUN.EQ.2>  CO  TO  C7 
1FIIFUN.E0.3I  CO  TO  08 
IF< IFIIN.E0.A1  GO  TO  69 
S3IJ)>FOCU$3-IFlC-FlXI«XNACin 
60  TO  17 

67  S1UI*F0CUS3-(F2C-F1X)*XNAG(II 
GO  TO  17 

66  S3U>«F0CUS3»IF1C«F1X)*XHA6<II 
GO  TO  17 

69  S)<J)«F0CUS1»(F*C«F1X>*XIUGH> 

17  V2< JI*(I3IJI-$21 J1IMFIASH3-FLASM2) 

V3  T  J  »-  <|3(J)-S1 ( J)>/<FLASH3-FLASH1I 
|F( 2X111. CO. 2X121)  GO  TO  70 
L3IJI-0 
GOTO  75 

70  131J1-1 
GOTO  80 

75  VOKJI-  VOlM)MFI*iHl/3.*(61**2«61*62*62**2)»FI*XH2/3.*IR2**2»62* 
1636039*2)) 

F16-F12 

h  VOL  1 I  .6  I  •  ri*XHl/3.*<61**2*61«02*62**2>*FI*XH2/3.*t62**2*62*63* 
163**2) 

XVOLIWJ)*  XVOLI lpj)*. 001*6 HO 

XL<6)>XI<6>»S06T((F1X'F2XI**2*1F1V-F2Y)**2)*XNA6II) 

OIAI ■M-DIAI  J)«2.*R2 
60  CONTINUE 

IMI.LT.NFLASH)  GOTO  169 
IF1NFLASH.E0.2)  GOTO  05 

1F(L1|J).E0.1.AN0.12I  JI.E0.1.AN0.L3I  D.EO.ll  GO  TO  169 
IFU116). EO.O. AND. L2IJ). EO.O.ANO. 1311). EO.O)  GO  TO  90 
IFI12IU.E0.1.AN0.L3I  D.EO.l)  GO  TO  95 
1FU1|J).E0.1.AN0.13IJ).E0.1)  GO  TO  95 
IF|llll).E0.1.AN0.L211).E0.1l  GO  TO  95 
TriL211I.EQ.0.AN0«L3l4).EQ*0)  GO  TO  100 
TF11U 2>.EO.O.ANO.L3IJ).EQ  .0)  GO  TO  100 
IFllll J).f 0.0. ANO. 1216). EO.O)  GO  TO  100 
65  lFUltll.E0.1.ANO.L211).E9.1)  GO  TO  169 
lFllltJI. EO.O. ANO. 1211). EO.O)  GO  TO  100 
IFIlllJI. EO.O. 06, 1212). EO.O)  GO  TO  95 
90  XNASSUHVOim*. 001/3. MHO 
XL  ID  I • XL  1 1 1  /  3. 

01AI3)«0tA(J)/3. 

GOTO  169 

95  XNASSI 6 ) •VOL IJI*.001*6HO 
GOTO  169 

100  XNAS$I2)«V0LIJI*. 001/2. «6H0 
XI 1 2 )*Xl (6 1/2. 

0IAI2 ) *01 A 12 )/2 • 

169  CONTINUE 
190  CONTINUE 


•612330 

•5 

6612330 

•  C 

•  612330 

•  7 

•612330 

•• 

•612330 

•9 

»6t2330 

90 

•612330 

91 

•RL2330 

92 

1612330 

93 

•612330 

9A 

•612330 

95 

•612330 

96 

•612330 

97 

•612330 

96 

•612330 

99 

•612330 

100 

•612330 

101 

6612330 

102 

•612330 

103 

•612330 

10* 

•612330 

105 

•612330 

106 

•612330 

107 

•612330 

lOt 

•612330 

109 

•612330 

110 

1612330 

111 

■612330 

112 

•612330 

113 

•612330 

M* 

S6L2330 

115 

•612330 

lie 

•RL2330 

117 

66L2330 

116 

B6L2330 

119 

■612330 

120 

R6L2330 

121 

*612330 

122 

*612330 

123 

•  612330 

12* 

6RL2330 

125 

•  612330 

126 

••12330 

127 

■612330 

126 

•612330 

129 

•612330 

130 

•612330 

131 

•612330 

132 

6612330 

133 

•612330 

13* 

6612330 

135 

6612330 

136 

ML  2330 

137 

••12330 

136 

•612330 

139 

•612330 

1*0 

••12330 

1*1 

■  612330 

1*2 

•612330 

1*3 

M12330 

1** 

•612330 

1*5 

••12330 

1*6 

••12330 

1*7 

ietmeak 

1 

60 


7 


c 

c 

c 

c 


8RGK  IS  THE  NUMBER  OF  JET  PARTICLES  THAT  HAVE  PARTICULATCO  A<  OF  JETBREAK 
EXPOSURE  TINES  1*  2.  AMO  3  RESPECTIVELY.  TDLAY  IS  THE  OELAY  TINE  JETPREAK 

BETWEEN  DETONATION  ANO  JET  FORNATION  INKROSEC).  JETBREAK 

JETBREAK 

REAOI5.PI  BROK1. BROK2.BRQK3. TOLAY  JETBREAK 

PZ-BROKl«BROKZ«SROK3  JETBREAK 

TZ»FIASH1pFLAJH2*FLA$H3  JETBREAK 

PBYT7«8*OKJ*FIA5H1PBROK2*FIASM2*BROK3*FLASMI  JETBREAK 

PBYP7.BROKl*»2.8ROK2P*2*RROKl**2  JETBREAK 

II.(PBYTI-IPIPTIT3.I  WIPBYFI-IPIPP2/3.II  JETBREAK 

PBAR*PZ/3.  JETBREAK 

TAVG«T2T3.  JETBREAK 

00  170  N»1»NPA*T  BRLE330 

IFINFIASH.ES. 2ISOTO  195  BRL2330 

VELINI*IVllNlPV2fN|PV3INM/3.P.l  BRL2330 

60T0  160  BRL2330 

155  VEL IN) •VIINI • , 1  BRL2330 

160  SUNL(N|.XL(N1»SUHL(N-1I  B>12330 

5UNNAS(N|.XHASS(NI»SUNHAS(N-1I  BRL2330 

ELOOINI-XllNI/OIAINI  BRLE330 

VELINI'VEKNIPIO.  BRL2330 

XKEINI-.9PXNA5S ( N ) P . 001* ( VE L (N ) P 1 000 . IPP2  BRL2330 

BURKE  INI >XKE IN) «SUHKE(N-1I  BRL2330 

PINI«VEL(N|PXNASS(NI  BRL2330 

TCTPINI-PINI  MOTP  (N-l  I  BRL2330 

SUNOI A  IN  I .  SURO I A  IN- 1 IpOIA(N)  RRL2330 

BREAK<NI«TAVS«f 1PIN-PBARI-T0LAY  JETBREAK 

GOTO  1  TO  BRL2330 

165  BREAK  111 *0.0  BRL2330 

1  TO  PRINT  lTl*  N.miNI.SUNLINI.BREAKINI  BRLE330 

1T1  F0*NAT(23X.12.1OX»F6.3»9X«F6.2»?X.F6.1l  BRL2330 

PRINT  92  BRL2330 

82  FORHAT I1H1.20X.' PARTICLE*. 6X> '  VEL0CITY1*»6X.  '  VE10CITY2'»6X»  BRLE330 

l  •  moc  I  m*. /,  22X. ''NUMBER*  »3X*»  I HH/NICROSEC I ’ .2* » • ( NN/NICROSEC I ' .  BRL2330 

22X.'»HH/HIC*05EC1'./I  9RL2330 

00  172  J.i.NPiRT  BRL2330 

172  PRINT  173.  J.V1IJI.V2IJI.V3UI  BRL2330 

173  F0RHAT<23X.T2.4X.F6.3I'"'»F6.3.9X»F6.11  BRL23JO 

PRINT  175  BRL2330 

175  FORMAT (1H1.20X. 'PARTICLE  LENGTH  OIA.  I/O  NASS  TOTAL  J  BRL2330 

1ET'./.22X. ’NURBER'.SX, • INNI’.AX/'  IHN!',9X.  MGR  ANSI',  2X,*  NASS  IGRANS  BRL2330 
2IN/J  BR  L  23  30 

DO  176  1*1 .HEART  BRL2330 

176  PRINT  ITT.  I.XLIII.DlMll.E;ODIII.TNASSIIItSUNNASIII  BRL2330 

177  F0RRA7(23Xvt2«RX.FA.lfAX.FA.1.2X.FA«I.AX.  F9.2.5X.F6.2J  BRL2330 

PRINT  603  BRL2330 

603  FORHAT (1H1.20X.  'PARTICLE '.*X« 'KASS1' >9X. • NASS2' .9X» • HASS3' «/» 22X.  BRL2330 

I  'NUNS  PR*  .AX.*  IGRANS  !'•  »  7X»*  ISRANB  J  •  »TX  ,  ■  I  GRANS!'  •  / 1  BRL2330 

DO  600  J'l.NPART  BRL2330 

1*1  8RL2330 

600  MR1TEI6.602I  I  J.  XVOl  1 1 .  J  I.  XVOl  II»1 .  Jl.XVOlt  1*2*  Jl  I  BRL2330 

602  FORHATI23X. 12. 6X. FB.6i6X.PS. A. 6X.F8.AI  BRL2330 

PRINT  1B0  BRL2330 

1B0  FORNATttHt.ZOX. 'PARTICLE  K.E.  TOTAL  JET  DISTANCE  FRON  CH  BRL2330 

URGE  BASE'  »  /  .  22X  »  'NUMBER  ' »  AX.  •  I  JOULES  I '  ,  3X»’  KE I JOULESI '  *2X.' FLASH  BRL2330 
21’. 21. 'FLASH  2'.2X.'FLASH  3«./l  BRL2330 

00  181  I«1.N*ART  SRI  2330 

181  PRINT  112.  UXKEIII.SUHKEIII.SIII  I .  S  2 1 1 1 .  S  31 1 1  BRL2330 

1*2  FORHATI23X.12. TXtFR.O.AX.FR.O.  AX.F5.0.  AX.  F9.0.AX.F9.0I  8RL2130 

PRINT  IBS  BRL2330 

183  F0RNATIIH1.20X. 'PARTICLE'. AX. ' NONENTUN'.AX, 'TOTAL  JET' . T.22X, 'NUMB  BRL2330 
IERP.AX.' IXG-N/SECI'. 3X.'N0NENTUN' *  J I  BRL2330 

00  189  W.NPART  BRL2330 

1*5  PR1NI  ISA.!. Pfll. TQTPtlt  8RL2330 


2 

3 

A 

5 

6 

7 

8 
9 

10 
11 
12 
13 
1  AB 
1A9 

150 

151 

152 
193 
15A 

155 

156 
1ST 
15* 

159 

160 
161 

IA 

162 

163 

16A 

165 

166 
167 
1 6* 

169 

170 

171 

172 

173 
17A 

175 

176 

177 
17* 

179 

180 
181 
1B2 
183 
1  BA 

185 

186 
1*7 
1 88 

189 

190 

191 

192 
143 
19A 

195 

196 

197 
19* 


61 


m 

F0RNATI21X*I2*9X*F6.Z»7X«F6.2) 

BRL2330 

199 

00  192  N»2*NPART 

BRL2330 

200 

SUNIENIN>«$UNLEN<N-1)«XLIN) 

BRL2310 

201 

X0IAIN>*X0IAIN-1>»0IAIN) 

BRL2330 

202 

IFIN.E0.2)  CO  TO  191 

BRL2330 

203 

DELVtN)«VELIN-ll-VEL<N) 

BRL2330 

209 

GO  TO  192 

BRL2330 

203 

191 

OELVtNI«VELtl)-VEllN) 

BRL2330 

206 

192 

SDELV1NI-SDEIVIN-1>*0EIVIN> 

BRL2330 

207 

00  193  J*3*NPART 

BRL23S0 

206 

195 

SUNDELtU>«DELVI*)»SUHDELU-l> 

SRI  2330 

209 

AVL  1*SUNHNPARTI  VFLOAT  INF  ART  1 

BRL2330 

210 

AVl2«SUNLENlNPART)/FlOATINPART-t) 

BRL23  30 

211 

AV01-SUH01AINPARTI/f LOATINPART) 

SRL2330 

212 

AV02-XDIAINPART > /FLOAT INPART-ll 

SRL2330 

213 

ADELV1*S0ELV INPART) /FLOAT  INPART-1 1 

BRL2330 

219 

A0ELV2*SUN0El<NP ART! /FLOAT <NFART-2I 

SRL2330 

213 

PRINT  200*  AVL1*  AVL2*  AV01*  AV02* AOELVlt ADELV2 

BRL2310 

216 

200 

FOX HAT  llHl#/7///*97X**WITH  JET  TIP  Nf 0  TT P*. /*20X»* AVER  AGE  PART 

BRL2330 

217 

1ICLE  L  ENCTH1 *7X*  F6.2*7X*F6.2*/*20X#VAVERAGE  PARTICLE  D!AHETER'»6X* 

BRL2330 

2  IB 

2F9.2*6X*F9.2*/»20X«*AVERAGE  CHANGE  IN  VEL0CITV*»6X*F9.2*9X*F9.2) 

SRL2310 

219 

PRINT  (09 

BRL2330 

220 

NRA*2 

SRL2330 

221 

NN*1 

BRL2330 

222 

TC*0 

BR12330 

223 

CALL  POLYlStVEl.Sl»NPART.AA,NRA,NN»CC*RR*AF,ERH$.SIC.TT*OET.IC I 

BRL2330 

229 

WRITEI6*91)  CCI1) 

BRL2330 

223 

osm-eem 

BRL2130 

226 

*1 

FORHAT 1 20X# •  VIRTUAL  ORIGIN  FOR  FLASH  1*<*F12.6) 

BRL2330 

227 

CALL  P0LVISIVEL»I2»NPART.AA.NRA,NN,CC.RR.AF,ERNS*SIG*TT»DET.  ICI 

BRL2310 

22B 

WRITE  1  6*  92)  CC(l) 

BRL2330 

229 

D3I2I-CCI1) 

BRL2330 

230 

♦  2 

F0RNATI20X,*  VIRTUAL  ORIGIN  FOR  FLASH  2*>*F12.6> 

BRL2330 

211 

IFIL.LT. 3)  GO  TO  915 

BRL2310 

212 

CALL  POLYLSI VEL*  S3«NPART«AA*NRA»NN*CCiRR*  AF*ERRS*SIG*TTtO£T*IC) 

BRL2130 

211 

VRITEI6#93)  CCI1) 

BRL2330 

239 

DSID’CCIl) 

BRL2330 

233 

91 

FOR HATI20X**  VIRTUAL  ORIGIN  FOR  FLASH  3«'»F12.6) 

BRL2330 

236 

c 

DRtFT 

l 

c 

ORIFT  CALCULATIONS 

DRIFT 

2 

c 

DRIFT 

1 

c 

ORIFT 

9 

c 

EACH  REF  REFERS  TO  THE  OISTANCE  BETWEEN  THE  REFERENCE  PARTICLE 

ORIFT 

3 

c 

AND  THE  FIOUCIAL  FRON  WHERE  ALL  OTHER  HEASURENENTS  ARE  TAKEN. 

DRTFT 

6 

c 

REF  tS  TAKEN  IN  HH  ANO  SHOULD  SE  THE  ACTUAL  DISTANCE  READ  OFF  OF 

ORIFT 

7 

c 

THE  FUN.  IT  WILL  THEN  SE  HULTIPLIEO  ST  THE  NAG  FACTORIWHtCH 

DRIFT 

6 

c 

SHOULD  HAVE  BEEN  ENTEREO  AS  DATA  IN  SRL2S30)  NOTE*  IF  XRA6  IS 

ORIFT 

9 

c 

LESS  THAN  1.  «  THEN  THE  FILH  PORTRAYS  AN  INAGE  LARGER  THAN  REALITY 

DRTFT 

10 

c 

DRIFT 

11 

REA0I5**)  REFI1)*REFI2).REFI3) 

DRIFT 

12 

REF  1 1 )  *REF  ID9XNAGI 1) 

ORIFT 

13 

REFI2)*REFI2)*XNAGI2> 

ORIFT 

19 

REFI3)<REF|3)*XNAGI3) 

ORIFT 

13 

DO  901  J*1  *3 

DRTFT 

16 

00  900  T*1*NPART 

DRIFT 

IT 

900 

A2U,I)«A2U*I)-REFU) 

ORIFT 

16 

901 

CONTINUE 

OR  TFT 

19 

PRINT  33 

ORIFT 

20 

33 

FORNATI 1H1  *21X*xPART  ICLE*  *  BX*  v  DEVIANCE  FRON  PATH  INHH*7«23X# 

DRIFT 

21 

l'NUNBER<»6X*<FLA3H  1  FLASH  2  FLASH  3*'/) 

DRIFT 

22 

DO  BOO  JJJ-l.NPART 

DRIFT 

21 

WRITE  16*36)  JJJ»AZtl*JJJI*AZI2*JJtfJ*AZI3*JJJ) 

DRIFT 

29 

36 

FORNATI 29X*I 2* BX* FP. 1*3X*F B.3*  3X»FB. 1) 

DRTFT 

29 

aoo 

CONTINUE 

DRIFT 

26 

00  491  1<1,NPART 

DRIFT 

27 

TiiiMSim-DsmifVELm 

ORIFT 

28 

T2II)*IS2III-0SI2))/VELLI> 

ORTFT 

29 

T3IT  >*IS3U)-DSI3))/Vflll> 

DRIFT 

30 

KD12*AZtl.I)*T2(I>/<AZI2,I  )*T1II) 1 

drift 

31 

K023.AZI2,1)*T3II)/C AZI3#t )*T2III  1 

DRIFT 

32 

K0ij-»m»n*Tjm/uz(3.i»*mti  i 

ORIFT 

33 

DO  409  8«1»3 

DRIFT 

34 

IFIA7IU#I)«REFIJ))  409 #402 #405 

ORIFT 

35 

402 

THETA(I>«0.0 

DRIFT 

34 

VETRAV(I)«0.0 

ORIFT 

37 

CO  TO  491 

DRIFT 

38 

409 

CONTINUE 

ORIFT 

39 

c 

DRIFT 

40 

c 

SOLVING  FOR  THETA  USING  FIRST  TVO  GOVERNING  EQUATIONS 

ORIFT 

41 

c 

ORIFT 

42 

TH12«ATANI2.0*K0127SQRT( 2. 01-1.0) 

ORIFT 

43 

IFITH12.lT.O.O)TH12«TH12*PT 

DRIFT 

44 

c 

ORTFT 

45 

c 

SOLVING  FOR  THETA  USING  LAST  TWO  GOVERNING  EQUATIONS 

ORTFT 

44 

c 

DRtFT 

47 

TH23#ATANI1.0-2.0/(SQRT(2.0)*K023n 

ORIFT 

48 

1F<TH23.LT.0.0)TH23-TH23*P1 

ORIFT 

49 

490 

IF IAIS ITH12-TH23 t.LT .81/2. 0)60  TO  449 

ORTFT 

50 

IFITH23.LT.TH12) GO  TO  448 

ORIFT 

51 

TH12#TH12#FT 

ORTFT 

52 

GO  TO  449 

ORIFT 

53 

448 

TH12«TH12-FI 

ORIFT 

54 

c 

DRIFT 

55 

c 

SOLVING  FOR  THETA  USING  FIRST  AND  LAST 

GOVERNING 

EQUATIONS 

DRIFT 

54 

c 

ORIFT 

57 

440 

T HI  3* A TAN I (KOI 3-1 .0)/f KOI 3 41.0) ) 

DRIFT 

98 

IF(TH13.1T.0.01TH13*TH13»PI 

DRTFT 

59 

470 

IF  I A8S (TH23-TH13 I.LT.FI/2.0)G0  TO  449 

ORIFT 

40 

IF  1 TH23.LT  «TH13)G0  TO  448 

ORIFT 

41 

TM3>TH13*P1 

ORIFT 

42 

CO  TO  449 

ORIFT 

43 

488 

TH13#TH13-PI 

DRIFT 

44 

c 

DRIFT 

45 

c 

FINDING  AN  AVERAGE  VALUE  NOR  THETA  ANO 

FLACING 

IT  IN  THE  PROFER 

DRIFT 

44 

c 

QUAORANT 

ORIFT 

47 

c 

ORTFT 

48 

449 

THV<TH12*TM23*TH13 1/3.0 

DRIFT 

49 

TTI I TH23-FI/2.0) 4AZI 2# I ) )  471,473,472 

ORTFT 

70 

471 

THETA! I 1-THV 

DRIFT 

71 

CO  TO  480 

ORIFT 

72 

472 

THETAI I I-THV4PT 

ORTFT 

73 

TH12#TH12#PI 

ORIFT 

74 

TH23«TH23»PI 

ORIFT 

75 

TH13#TH134F1 

DRIFT 

74 

IF  I THETAI I )-THV. GT .1 .7)THV#THV4FI 

ORIFT 

77 

GO  TO  480 

DRIFT 

78 

471 

TFITH13)  474, 479# 475 

ORIFT 

79 

474 

THETAI I )*THV#FI 

DRIFT 

80 

TH12-TH124FI 

ORIFT 

81 

TH23>TH23»PI 

DRTFT 

82 

TH13-TH134FI 

ORIFT 

83 

IF ITHETAI 1 ) -THV . GT . 1 .7) TNV*THV*PI 

DRIFT 

84 

GO  TO  480 

DRTFT 

85 

479 

THETAI  D.THV 

DRIFT 

84 

c 

ORIFT 

87 

c 

CALCULATING  RELATIVE  TRANSVERSE  VELOCITIES  UN  HE 

ASSIGNING 

ORIFT 

88 

c 

VE1CHTE0  VALUES... IOEFENOING  ON  THE  FARTICLES* 

ORIENTATION 

DRIFT 

89 

c 

V«TH  THE  FLASH  R-RAV) 

ORIFT 

90 

<,80  V£TR12.AZ(l,tlMCOS(TNI2-PI/4.0l*Tl<!MR1000.0 
VE  TR23.AZI2.  I  IflC0SITH23)«  Tim  1*1000.0 
VETR13-AZ ( 3.  !»/ IC0S<TM13.PI/4.0I*T Sill  1*1000.0 
E  PJ(1.  I  I*»ZI  1«I  >*SlN(TM12-H/4. 0IR(TH12-TMVI/Tllt  1*1000.0 
I  /IC0S(TH12-PI/*.0>>*»2 

fP3I2.1>*AZI2.I)«SIN(TH23)*(TH23-TMVl*1000.0/T2tI> 

1  TIC  0S(  TH23  I  I  **2 

EP$tI»ll«AZ«l»II*SIN(TH13*PI/4.0I*ITH13-TMV>/T3III*l  000.0 
1  7(C0SITH13+PI/4.OII«*2 

FP5I4,1)-AZ<  1. II *SIN( TH13-PI /4.0 I* (TM13-THV I / 11(11*1000.0 
1  MC0S(TH13-ft/4.0ll**2 

EPS(5. 1 l»AZ 1 2.1) *SIN(TH12)*(TM12-TMV 1*1000.0/72(1 1 
1  MC0S(TH12II**2 

EPS(6»n*AZI3.D*SIN(TH23*RS/4.01*(TH23-THVI/T3(f  1*1000.0 
1  /( C0S(TH23.PI/4.0)I**2 

DO  481  J-1.6 

481  EPS(J#I>**BS(EPS(J'»1)> 

EPSUH«0.0 

DO  482  *»1.6 

482  EPSUH.1.0/EPS(J»I)*EFSUN 
DO  483  J»l»6 

483  WT(JI*1.0/(EPS(»»I)*EPSUNI 

VETRAV(TI*NTI1I*VETR12.VT(2I*VETR23*WT(3I«VETR114WT(4>*VETR12 
1  ♦VT(*1*VETR23.VTI61*V£TR11 

CONVERT  THETA  TO  DEGREES 

490  THETA(II-180.0*TNETA(I)/PI 

491  CONTINUE 
PR1HT  ^92 

492  F0RNAT(1H1.,21X,*RARTICIE'.9X, 'RELATIVE  TRANS VERSE*»5X, 'RELATIVE  • 

1  *  DEVIANCE  *-./. 23X*'NUN  BE  R*.8X»*  VELOCITY  IN/SI  *.9X.  •ANGIE*. 

>  •  (OEGREES)1. /I 

DO  493  K-l.NPART 

WRITE  (6.  494 1  (C,VEr*AV(Kl.r«ETA(*I 

494  P0RHATIZ4*,I2.1»X>F6.2.17X.FS. 1) 

493  CONTINUE 

TRANSLATE  APPARENT  VALUES  INTO  ABSOLUTE  VALUES 


GUESSING  AT  AN  INITIAL  OEVIANCE  ANGLE  FOR  THE  REFERENCE 

TNIN>THETA(U 

TNAX-THETAU) 

HINOEXM 

NAXOEX.l 

00  391  W.NPART 

lE(THETAm.LT.TNIN.0*.TNlH-THETA(I).lT.-300.1S0  TO  392 
GO  TO  393 

392  THIN«TMETA(  I 1 
NINOEX.I 

GO  TO  391 

393  IF(TMETAUI.GT.TRAX. OR. THAX-THETAdl.GT. 300.1  GO  TO  394 
CO  TO  391 

394  TNAX>THETA(I) 

NAXDEX-I 

391  CONTINUE 

IMNINDEX.CT.NAXOEXI  CO  TO  399 

1NCR..003490T 

THETA0-THETKHAX0EXI-180. 

GO  TO  396 


OR  TFT 

91 

ORIET 

92 

DRIFT 

93 

ORIFT 

94 

DRIFT 

95 

DRIFT 

96 

DRIFT 

97 

OR  TFT 

98 

DRIFT 

99 

OXtFT 

100 

ORIFT 

101 

DRIFT 

102 

DRIFT 

1C3 

ORIFT 

104 

ORIFT 

103 

DRIFT 

106 

DRIFT 

107 

DRIFT 

108 

ORIFT 

109 

DRIFT 

110 

DRIFT 

111 

DRIFT 

112 

DRIFT 

113 

drift 

114 

drift 

115 

ORIFT 

116 

DRIFT 

117 

ORIFT 

118 

DRIFT 

119 

DRIFT 

120 

ORIFT 

121 

DRIFT 

122 

ORIFT 

123 

ORIFT 

124 

ORIFT 

125 

DRIFT 

126 

ORIFT 

127 

ORIFT 

128 

DRIFT 

129 

ORIFT 

130 

ORIFT 

131 

ORIFT 

132 

ORIFT 

133 

DRIFT 

134 

ORIFT 

135 

ORIFT 

136 

ORIFT 

137 

DRIFT 

138 

ORIFT 

139 

ORIFT 

140 

DRIFT 

141 

ORIFT 

142 

DRIFT 

143 

ORIFT 

144 

ORIFT 

145 

ORIFT 

146 

ORIFT 

147 

ORIFT 

148 

ORIFT 

149 

ORIFT 

150 

ORIFT 

1*1 

ORIFT 

152 

DRIFT 

153 

ORIFT 

154 

395 

INCH— .0039907 

DRIFT 

155 

THETA0-THETAIHINDEXI-180. 

DRTFT 

136 

396 

ORCS-O.O 

DRIFT 

1ST 

00  399  I-l.NFART 

DRIFT 

1  SB 

THETA 1 1  l-TMETAM  >9PI/ 180.0 

DRIFT 

139 

VXI I )-VETRAV 1 I )*C0I ITHETAI  111 

OR  I  FT 

160 

399 

VVlII-VETRAVUMSINITHETAdll 

DRIFT 

161 

THO»-THETAO*PI/1»0.0 

DRIFT 

162 

902 

llll'llil 

ORTFT 

163 

ITOTL-IF IX 1 • 122173/3(5 (IOC* • 1 

DRIFT 

169 

00  1100  NOT-l-ITOTl 

DRTFT 

163 

T HOP - THOR* INCR 

ORIFT 

166 

903 

VTO-O.O 

DRIFT 

167 

t 

DRIFT 

16B 

c 

FINDING  THE  HACN1TU0E  OF  REFERENCE  TRANSVERSE  VEIOCITTTHAT  MILL 

DRIFT 

169 

c 

HAKE  IT  FROFORTIONAl  TO  THE  AXIAL  VELOCITV 

DRIFT 

170 

t 

DRIFT 

171 

00  736  I-l-NPART 

ORIFT 

172 

736 

VTABSII l-VETRAVII 1 

DRIFT 

173 

909 

VTO-VTO*10.0 

DRIFT 

179 

INC-100 

DRIFT 

173 

GO  TO  909 

DRIFT 

176 

905 

VTO-VTOvl.O 

DRIFT 

177 

INC-10 

DRIFT 

170 

GO  TO  909 

ORIFT 

179 

906 

VTO-VTO-.l 

DRIFT 

1  BO 

INC-1 

ORIFT 

1 01 

909 

VXO-VTOFCOSI THOR  1 

ORIFT 

102 

V70-VT09SINITH0R) 

ORIFT 

103 

DO  907  l-l-NPART 

DRIFT 

1  09 

907 

VTABSII) -SORT  1 1 VI  II)-VX0>«*2MVTMI-V»0t*9«l 

ORIFT 

1B3 

c 

DRIFT 

106 

c 

STATISTICAL  calculations  FOR  A 

LINfAR  REGRESSION 

DRIFT 

187 

c 

DRIFT 

IBB 

002 

TS-0.0 

DRIFT 

189 

AS-0.0 

DRIFT 

190 

TBVAS-0.0 

ORIFT 

191 

TRVTS-0.0 

ORIFT 

192 

ARY AS-0.0 

DRIFT 

193 

KlESS-0 

DRIFT 

199 

DO  900  1-t.NPART 

ORIFT 

195 

IFI ARSI VTARSTI l-VTOI  .IE.  1.01 

CO  TO  700 

DRIFT 

196 

TS-VTARSIII-TS 

DRIFT 

197 

AS-VEUEI*AS 

ORIFT 

1 90 

TRVAS-VTARSII l*VELII 1 *TBVAS 

DRIFT 

199 

T8TTS* VTABSII 1*9 2*T3TTS 

ORIFT 

200 

ABVAS-VEL 1 I )9*29 ABYAS 

ORIFT 

201 

CO  TO  900 

ORIFT 

202 

700 

KLESS-KLESS*! 

DRIFT 

203 

900 

CONTINUE 

ORIFT 

209 

C01- CTBT AS- (TS»A*/ (NPART-XLESS 1 II / ITBTTS- ITS992/ 1 NFART-KIESSI 1 1 

ORIFT 

205 

C02-ITRVAS-I TS9AS7INFART-NLESSI 1 1  /  I  AIVA5-IAS992MNPART-KLESSI 1 1 

DRIFT 

206 

C0-IC01-1. 07 C02 172.0 

ORIFT 

207 

TBAR-TS/INPART-KLFSS) 

DRIFT 

208 

ABAR- ASM  NFART-KLESS  1 

ORIFT 

209 

X1NT-TRAR— ARAR/CO 

DRIFT 

210 

I9IXINT.lt. 0.0  .AND.  INC.E0.100)  GO  TO  909 

DRIFT 

211 

091X1 NT. LT.0.0  .AND.  INC. EQ. 101 

60  TO  905 

ORIFT 

212 

IFIXINT.lT.0.0  .ANO.  INC.EO.ll 

CO  TO  906 

DRIFT 

213 

EFIINC-10)  012-911-910 

ORIFT 

219 

910 

VTO-VT 0-11.0 

ORIFT 

213 

GO  TO  90S 

DRIFT 

216 

911 

VTO-VTO-1 . 1 

ORIFT 

217 

GO  TO  906 

ORIFT 

2  IB 

65 


or»or*r>r»r*oo 


912 

RI2)«SORTIC01*C02) 

DRIFT 

219 

lFIRI2l.tT.RUII  tO  TO  1101 

DRIFT 

2  ZO 

tO  TO  1100 

DRIFT 

2  21 

1101 

Rl 1 )*R 12) 

DRIFT 

222 

VTOB-VTO 

DRIFT 

2  23 

THOR »• THOR 

ORIFT 

229 

00  1'102  II«1,NFART 

DRIFT 

223 

1102 

VTIII6),VTABSI16) 

ORIFT 

226 

1100 

CONTINUE 

DRIFT 

227 

VTO«VTOB 

DRIFT 

226 

TMOR-THORB 

ORIFT 

229 

DO  1103  181*1, NFART 

DRIFT 

230 

1103 

VTABSIIBII-VTBUII) 

ORIFT 

2  31 

190 

DO  (13  N*1»NFART 

ORIFT 

2  32 

IFI I VTABSIN l-VTO 3 . L E .  1.0)  GO  TO  811 

ORIFT 

233 

FH1IN)*THETAIN)-TH0R 

DRIFT 

239 

BTTAIN>*ASINIVETRAVIN)*SIN{FHIIN>)/VTAtStNI)*180.0/FI 

DRIFT 

233 

A VS TIN)* VE11N) /VTABSIN 1*1000.0 

ORIFT 

236 

tO  TO  813 

DRIFT 

237 

Oil 

VTABSIN>*0.0 

ORIFT 

238 

BETA(N)*0.0 

DRIFT 

2  39 

AVSTINI-0.0 

DRIFT 

290 

913 

CONTINUE 

DRIFT 

291 

FR1NT  702 

ORIFT 

292 

702 

FORNAT 1 1H1 , 21X FART1CLE  *, 3X,*  ABSOLUTE  TRANS  VERSE*. 3X. • ABSOLUTE  *, 

DRIFT 

293 

1 

1  m ANGULAR * , 3X«  *  AXI AL  VELOCITY  OVER* ,/. 23X, 'NUHBER»,9X»* VEL ', 

DRIFT 

299 

2  •0C1TY  IM/S>',9X, 'DEVIANCE  WRT  REFERENCE  TRANSVERSE  VELOC', 

DRIFT 

295 

3  MTT't/1 

DRIFT 

296 

DO  710  J*l, NFART 

DRIFT 

297 

WRITE (6* 7071  J, VTABS II l>8E TAI J I, AVSTI J I 

ORIFT 

298 

707 

F0RHATI29I,I2,16X,F6.2,13S,F6.1,1TX,F6.2> 

ORIFT 

299 

710 

CONTINUE 

DRIFT 

230 

FRINT  712, VTO 

DRIFT 

231 

712 

FORMAT  If, 221, •THE  TRANSVERSE  VELOCITY  OF  THE  REFERENCE  IS  '» 

ORIFT 

232 

1  F6.2-,’  N/SU 

ORIFT 

233 

FRINT  713,  RI2I 

ORIFT 

239 

713 

DRMAT I /,22X, *  THE  CORRELATION  COEFFICIENT  FOR  THIS  REGRESSION*, 

ORIFT 

233 

1  •  IS  *,F3.9) 

DRIFT 

236 

C 

FENJET 

1 

c 

FENJET 

2 

c 

FENJET 

3 

c 

FENJET 

4 

c 

ECONST  IS  THE  ENERGY  REQUIRED  TO  CREATE  ONE  CUBIC-NH  OF  HOLE 

FENJET 

5 

c 

VOLUME.  RHOJ  It  JET  DENSITY  IN  t/CU-HN.  SIGJ  IS  JET  STRENfTH  IN 

FENJET 

6 

c 

MEGAFASCALS..  RHOT  AND  SIGT  ARE  THE  COBRESFONDTNG  TARGET 

FENJET 

7 

c 

VARIABLES 

FENJET 

8 

c 

FENJET 

9 

READI3, 300 )ECONST, RHOJ, SIGJ, RHOT, SIGT 

FENJET 

10 

300 

FORMAT I3F10. A ) 

FENJET 

11 

c 

FENJET 

12 

c 

THIN  IS  THE  MINIMUM  STAHOOFF  AT  WHICH  FENJET  HILL  EVALUATE  THE 

FENJET 

13 

c 

GIVEN  DRIFT  DISTRIBUTION.  IHAX  IS  THE  MAXIMUM  STANOOFF  TO  BE 

FENJET 

19 

c 

ANALYZED.  ZINC  IS  THE  INCREMENT  BETWEEN  STANDOFF  ANAL7ZATI0N. 

FENJET 

13 

c 

All  MEASUREMENTS  ARE  IN  NH. 

FENJET 

16 

c 

FENJET 

17 

READI3, 3971 ZMIN.ZMAX, ZINC 

FENJET 

18 

397 

F0FMATI3F10.il 

FENJET 

19 

DO  301  I*1»NFART 

FENJET 

20 

BETA!  1 )*6ETA<1>  *3.19139/180. 

FENJET 

21 

RAOIUS 1 I 1 *0. 0 

FENJET 

22 

FENI 1 1 *0.0 

FENJET 

23 

Rill), XI II) 

FENJt  r 

29 

RVEll I )*VELI I ) 

FENJET 

23 

301 

CONTINUE 

FENJET 

26 

66 


NNl*NPART-l 

PENJET 

27 

GANNA* SORT IRHOT/RHOJ 1 

PEN1ET 

21 

2SU60*ZHIN-ZIHC 

PENJET 

20 

K2SUH*0 

PENJET 

30 

260 

ISU60*Z9Ut0»ZINC 

PENJET 

31 

KZSUH*K2SUH»1 

PENJET 

32 

P£NNAXIKZSUNt«0. 

PENJET 

33 

ZSUBIKZSUNI*ZSUBO 

PENJET 

36 

c 

PENJET 

3* 

c 

THIS  1007  EVALUATE J  EACH  0IC1TI2ED  PAST JCIB 

PENJET 

36 

c 

PENJET 

37 

201 

00  200  f*l.NPART 

PENJET 

3B 

NCCII1*0 

PENJET 

39 

N0»0 

PENJET 

AO 

260 

TI 1 1* I  ZSUBO-Xl  I !  l-IDS (1I-D3I 21 »0S 13) 1/3. -2.99*1 7VELI  It 

PENJET 

AI 

0ELT*3./VELtll 

PENJET 

A2 

JNIN-l 

PENJET 

A3 

c 

PENJET 

AA 

c 

LOCATING  A  PARTICLE  IN  3-0  SPACE. ..THEN  FINDING  WHETHER  IT 

STRIKES 

PENJET 

A5 

c 

THE  HOLE  PROFILE.  IF  IT  DOESN’T^  THEN  TINE  IS  INCRENENTED 

UNTIL 

PENJET 

A6 

c 

INPACT  OCCURS. 

PENJET 

A7 

c 

PENJET 

AS 

201 

T  ( I l*T (1 1 »OElT 

PENJET 

A9 

20* 

XI  l)*mHl/AVST  <1 19T<1IPC0S<BETA<III 

PENJET 

*0 

Y<II*K<E|9TAN<BETA<II1 

PENJET 

*1 

Z<I  l*YEL<I)9T< I  l-ZSU60*KL<  II* <0E<II<0S<ZM0S<3I< /I. 

PENJET 

*2 

IFIN0.E0.2t  GO  TO  210 

PENJET 

*3 

IFIZI1I.LT. 0.01  60  TO  201 

PENJET 

SA 

1FK.E0.1  .AND.  Till. LT. BREAK (1 1 1  GOTO  202 

PENJET 

** 

IFII.EO.lt GO  TO  210 

PENJET 

96 

KLESS*I-1 

PENJET 

*7 

00  206  8* JNIN.KL  ESS 

PENJET 

*8 

IF(*ENIII.IT.. 01160  TO  206 

PENJET 

*9 

V<I«lt*X<It-X<Jl 

PENJET 

60 

»(I.2fT(ll-TUI 

PENJET 

61 

V<1>3I*ZM  l-ZMI 

PENJET 

62 

VNI It- SORT  IP  1 Irl t**2*» II, 2 t **2*¥<I.1IP*2I 

PENJET 

63 

VAN<1l*V<l.ll*GRA0<J.II*V<I.2l*GRA0<J'ZI*VtI»3|9GRAD<J»ll 

PENJET 

6A 

VANlIfABSIYANlIt) 

PENJET 

6* 

VPN(TI«IORT(ABS(VN(IIP*2-VAHIIt**2tl 

PENJET 

66 

IF|  IVPNI  It  *01  AM  )/2. 1»*2/RA0IUSI  J  )**2*VAN<II  »P2/PEN(  J)*»2. 

GT.l.OI 

PENJET 

67 

1 

I  GO  TO  206 

PENJET 

6B 

IFIN0.E0.1I60  TO  210 

PENJET 

69 

JNIN*8 

PENJET 

70 

GO  TO  201 

PENJET 

71 

206 

CONTINUE 

PENJET 

72 

207 

N0*1 

PENJET 

73 

TFIJNTN.EO.KLESS  .AND.  ZIII.LE..1I  GO  TO  206 

PENJET 

7A 

T<I1*TII t-OELT/6.0 

PENJET 

75 

GO  TO  20* 

PENJET 

76 

200 

T<ll*mUSO-XL<II-<OS<ll*OS<2l*OS<3l  t/3.1  /VELUt 

PENJET 

77 

NO-2 

PENJET 

78 

GO  TO  20* 

PENJET 

79 

c 

PENJET 

80 

c 

FINDING  THE  EFFECTS  OF  THE  PARTICLE  PENETRATION 

PENJET 

81 

c 

PENJET 

82 

210 

IFII.NE.l  .AND.  NO.NE .2 t GO  TO  202 

PENJET 

83 

GRA0I1* 11*0.0 

PENJET 

8* 

CRAOU. 21*0.0 

PENJET 

8* 

6RA0<1. 31*1.0 

PENJET 

86 

GO  TO  261 

PENJET 

87 

c 

PENJET 

88 

c 

IF  A  PARTICLE  TS  NOT  IN  THE  CONTINUOUS  NODI.  THE  AXIS  OF 

PENJET 

89 

c 

PENETRATION  IS  EVALUATEO 

PENJET 

90 

67 


c 

penjet 

91 

202 

TF(T(tl.LT.8REAK(II  INCC(I)  >1 

PENJET 

92 

IFII.E0.1  .AND.  NCCd  I.EQ.  1160TC1999 

penjet 

93 

IFtT(1l.lT.8REAK(I))Xl(  I)*Xl(l)*Tdl/8REAI<d) 

PENJET 

99 

im*i-i 

PENJET 

93 

IF  (T(I ). IT. BREAK  U)  .AND.  211).  ST. 2(1  XI)  (GOTO  1999 

penjet 

96 

00  209  LJ*l#3 

PENJET 

97 

VAU,LJI*6RA0(JHIN»LJI*VANdl 

PENJET 

98 

VFC  t.lJI*VII*LJI-VAI!fl8) 

PENJET 

99 

AR<l,LJ>*VA(IfLJ)/VAN(II 

PENJET 

100 

209 

PRII»IJ)*VPII»IJI/VPN(II 

PENJET 

101 

RGRA0(I*1I-VPN(II*2.7RA0IUS(JHIN)**2 

PENJET 

102 

PGRAOU. 21-0.0 

PENJET 

103 

RGRA0(l»3>-VAHd  l*2./PEN(  JN1NI*«2 

PENJET 

109 

R6N(II-SORT(PGRAO(I» I 1 »»2* RGRAD (I » 31 *»2 1 

PENJET 

103 

RGRADI  1.1 1 -RGRADII.il/RGNIl  I 

PENJET 

106 

RGRAD(1»3)-RGRAD(  U3)  /RGH(I) 

PFNJET 

107 

DO  219  LJ-1,3 

PENJET 

108 

219 

GRAD!  I#LJ)-RGRAD(I»1>*PR1I.IJ1«RGRA0(1*3I*AR(I»IJ) 

PENJET 

109 

GN-50RTIGRADII. 1)**2«GRA0( I- 2 ) *»2«GR AO ( 1. 3 )♦ *2 ) 

PFNJET 

110 

00  219  U«l,3 

PENJET 

111 

219 

GRAOINl  JI*GRAO(l.LJ)/GH 

PENJET 

11? 

GOTO  281 

PENJET 

113 

c 

PENJET 

119 

c 

CONTINUOUS  PENETRATION  AXIS  U  A10N6  P ARTICLE  FLIGHT  ARIS. 

PENJET 

113 

c 

PENJET 

116 

19*9 

VECT"SQRT(X(ll**2*T(l)**2*(2(I)*ZSUB0-(0S(l)*DS(2l*DS(3l)/3-)**2) 

PENJET 

117 

GRAOII>ll*X(t)/VECT 

PENJET 

118 

6RAD(I.2)*Y(]I /VTCT 

PENJET 

119 

GRA0II-.31-U  U)*2SU6O-(0Sm*DS<2l*0Sml/3.)/VECT 

PENJET 

120 

281 

PENNAX (K28UN )*XLd)/GANNA*PENNAX(KZSUHI 

PENJET 

121 

1F( I.EO-1  .OR.  NO. EO. 2  .OR.  NCC ( I > .EO.l IGOTO  262 

PENJET 

122 

C 

PENJET 

123 

C 

EFFECTIVE  LENGTH  AND  EFFECTIVE  VELOCITT  CALCULATIONS. 

PENJET 

129 

C 

PENJET 

123 

XL(T)-GRAD(I.3)*XL(1> *SORT (GRAD(lfl)**2*GRA0(If2)**2) 

PENJET 

126 

1  *01911) 

PENJET 

127 

VEL 1 1 >• GRAO ( I »  3)«VEL (II 

PENJET 

128 

c 

PENJET 

129 

c 

PENETRATION  CALCULATIONS  AS  A  FUNCTION  OF  EFFECTIVE  VELOCITY. 

PENJET 

130 

c 

PENJET 

131 

IF ( VEL ( I )• GT •  .8)  GO  TO  260 

PENJET 

132 

PENdl-O.O 

PENJET 

133 

RADIUS ( I 1*0.0 

PENJET 

139 

Udl-0.0 

PENJET 

133 

HVOLd  1*0.0 

penjet 

136 

GO  TO  209 

PENJET 

137 

280 

1F(VEl(II.GT.1.9l  60  TO  261 

PENJET 

138 

U(I>*0.0 

PENJET 

139 

HVOLd  )*XKE(1  l/ECONST 

PENJET 

190 

PEN(II*(VEl(ll-.ai**2/.72«Xldl 

PENJET 

191 

RADIUS (I )*SORTI 3. *HVOL ( I )  / (2. *3.191 39*PEN(1 1 1 ) 

PENJET 

192 

GO  TO  209 

PENJET 

193 

201 

1F(VEL((I.GT.2.)G0  TO  262 

PENJET 

199 

U(II*0.0 

PENJET 

193 

HVOL( l l*XKE ( I l/ECONST 

PENJET 

196 

PEN! II* (-(VEL (I 1-2. )**2/.72*l,l*Xl(I 1 

PENJET 

197 

RA0TUSdt*S0Rr<3.*HV0l(II  /(2.*3.191S9*PENdlll 

PENJET 

198 

GO  TO  209 

PENJET 

199 

282 

RAD-30RT(VELdl**2*RH0J*RH0T-2.*(RH0J-RHQT)*(51GJ-SIGT)/l.E6) 

PENJET 

150 

Ud  1*1 RPOJ*VEl(( l-RAO  l  /( RH0J-RHCT1 

PENJET 

131 

PEN(I I *(RHOJ*VEL (I)*Xl(I)-XL(I) *RADI 2( RA0-RH0T*VEl( III 

PENJET 

132 

HVOL (1 >*XKE( I l/ECONST 

PENJET 

133 

RAD(US(1I*SQRT(  3,*HV0L  ( 1 1  /  (2  .*!.  19l59*PENd  1 1 1 

PFNJET 

139 

209  CONTINUE 

EENJET 

133 

279  WRITE | 6,211) 

PENJET 

136 

211  FORMAT  1 1H1.10X. ' PARTICLE ' • 5X*' HOLE '»6X, • IMPACT 1 , SI, ' PENETRATION  *, 

PENJET 

157 

l  'PENETRATION', 6X, 'STATUS '»<P12X» 'NUMBER ',3  X,*  VOLUME'.  6X»  *TIMEE#17X» 

PEN  JET 

138 

2  •VEtOCITT,»/,2J*«  MCU.MM)  1HICROSEC  )•  ,  ♦*,  •  t  MM)  SX.'  INHF'  , 

PENJET 

139 

3  •MICROSEC)',/) 

PENJET 

160 

00  711  1*1 ,  NPAR T 

PENJET 

161 

TFINCCUI.EO.ltGOTO  281 

PENJET 

162 

NR  ITEI 6, 212)1  .HVOLII  )*TI  D.PENII  )»Ult  ),MPA,MRAR,MRART 

PENJET 

163 

GO  TO  213 

PENJET 

164 

281  WRITE  (6,  212)  WHVOl 1 1  >•  Tl II*PEN  11 1.UITI , HCQ,NCON,NCONT 

PENJET 

165 

212  F0RHATU4K,I2,7X,F6.0,4X,FS.l,7X,F4.neX,FS.SR7Xv8A4l 

PENJET 

166 

213  CONTINUE 

EENJET 

167 

EMAX-0.0 

PENJET 

168 

00  213  1>1,MPART 

PENJET 

169 

IF(ZU).LF.PNAX)  GO  TO  213 

PENJET 

170 

PMAX.2U  1 

PENJET 

171 

0PNAX.PENU)*GRADII,3) 

PENJET 

172 

213  CONTINUE 

PENJET 

173 

PEN3UN>PNAX,0PHAX 

PENJET 

174 

PENT IKZ5UH). PENS UN 

PENJET 

173 

WR1TE(6,216)PEN]UN 

PENJET 

176 

216  F0RMATI1N0,10X,'THE  TOTAL  PENETRATION  PREDICTED  8Y  THIS  •» 

PENJET 

177 

1  'PROGRAM  IS  »»F6.1,«  NNM 

PENJET 

178 

PRINT  237,  ZSU80 

EENJET 

179 

237  FORMAT  I1H0»10X, 'THESE  VALUES  WERE  CALCULATED  AT  A  STANDOFF  ', 

PENJET 

180 

1  'OF  ' ,  F6.0,  •  NM.') 

PENJET 

181 

VRTTEI6.220) 

PENJET 

182 

220  FORMAT  UM1 ,14X,«  X, T ,  AND  t  COORDINATES  PENETRATION  NAX  CRA', 

PENJET 

183 

1  *  TER  RADIUS  AXIAL  DIRECTION  OF  PENETRATION') 

PENJET 

184 

00  218  1>1,NPART 

PENJET 

183 

NRl  TEI6,  21 711, XIII.Yttl.ZIl), PENH ),  RADIUS  It  ),GRAD(  1,1  )  , 

EENJET 

186 

1  6R AO  11,2), GRAD It, 3) 

EENJET 

187 

217  F0RNArmx,I2,3F10.2,3F10.4l 

EENJET 

188 

218  CONTINUE 

PENJET 

189 

DO  263  1-1, HEART 

EENJET 

190 

VEL II) .RVELIt) 

PENJET 

191 

KWAl  It  1*0 

NJET 

192 

263  XLI I ).RL II ) 

PENJET 

193 

IEI2SUB0.lt. IMAX)  GO  TO  230 

PENJET 

194 

WRITEI6.290I 

PENJET 

193 

290  F QRMATI 1H1 , '  VIRTUAL  ORIGIN  ZOINHI  PENETRATION! HR)  MAX  ', 

PENJET 

196 

♦  'PENETRATION INN) •*/) 

PENJET 

197 

DO  289  l.l.KZSUN 

EENJET 

198 

289  WRITE  1 6, 292) ZSUR III, PENT  11 I.PENMAXII) 

PENJET 

199 

292  F0PMATI10X,FS.O,T32,F7.2,T33«F7.ZI 

PENJET 

200 

313  CONTINUE 

PENJET 

201 

303  FORMAT  1 1H1 ) 

8PL  2330 

2  37 

00  310  JN>  1, ME ART 

BR12330 

2  38 

VOLUNI-O. 

8R12330 

2  39 

XKJNI-O. 

BEL  2330 

240 

310  OTAIJNI.O. 

6RL2330 

241 

300  CONTINUE 

BRL2330 

242 

STOP 

8R12330 

243 

ENO 

BRL2330 

2  44 

69 


"JETBREAK" 


BREAK (N)  -  breakup  time  of  particle  N  in  psec 

1) 

BROK  2  >  -  observed  number  of  jet  particles  already  particulated  as  of 
3)  1 

exposure  2 
3 


l)  1 

FLASH  2>  -  time  after  detonation  of  round  for  exposure  of  film  2 
3)  3 


PBAR 

PBYPZ 

PBYTZ 

TAVG 


statistical  variables  used  in  calculation  of  breakup  time 


TDLAY  -  time  delay  between  detonation  of  round  and  formation  of  jet 
Z1  .  -  time  between  breakup  of  adjacent  jet  particles 


Figure  Al.  A  Diagrammatic  Sketch  Identifying  Some  DRIFT  Variables 


"DRIFT" 


AVST  (I) 


AZ  (1,  I)) 
AZ  (2,  I)J 
AZ  (3,  I)) 


BETA  (I) 

DS  (I)) 

DS  (2)  S 
DS  (3)) 

EPS  (1,  I) 
EPS  (2,  I) 
EPS  (3,  I) 
EPS  (4,  I) 
EPS  (5,  I) 


EPSUM 


KD  12  1 
KD  23  > 
KD  13  ) 

MAXDEX 

MINDEX 

PHI  (I) 


R  (2)1 

R  (Dl 


51  (1)1 

52  (I)} 

53  (I)) 


-  the  axial  velocity  divided  by  the  transverse  velocity  for 
particle  I 


-  deviance  of  particle  I  with  respect  to  the  reference  on 

1 

exposure  2 
3 

-  absolute  angle  formed  between  drifting  reference  and  particle  I 

1 

-  virtual  origin  calculations  based  on  exposure  2 

3 


-  estimate  for  sensitivity  of  transverse  velocity  V  with  respect  to 

dV 

small  changes  in  drift  angle  8  (of  the  form  A9)  using  all 

combinations  of  drift  velocities  and  drift  angles  calculated  from 
the  3  radiographs 


l  eric  Tv — tT  a  constant  used  to  scale  the  EPS  values  to  their 
K=1  *■  ’  ’  proper  magnitude 

1  2 

-  a  constant  relating  time  delay  between  exposures  2  and  3  used  in 

drift  velocity  calculations  1  3 

-  the  particle  #  that  satisfies  TMAX 

-  the  particle  #  that  satisfies  TMIN 

-  the  supplementary  of  the  relative  drift  angle  between  particle  I 
and  the  reference 

-  correlation  coefficient  for  V.  vs.  V_  plot  for  reference  drift 
angle  presently  being  analyzed. 

previously 


1 

-  axial  distance  from  charge  base  to  particle  I  on  exposure  2 

3 


T1  (I)) 
T2  (I)> 
T3  (I)) 


THETAO 


TMAX 

TWIN 

TH12) 
TH23  > 
TH13) 

THV 

THETA  (D 


THOR 

VETR12) 
VETR23  > 
VETR13) 

VETRAV  (D 


VTABS  (D 

VTO 

VXO 

VYO 

VX  CD 
VY  CD 

WT  (I) 


ABAR  * 

ABYAS 

AS 

CO 

COl 

C02 

TBAR 

TBYAS 

TBYTS 

TS 


\ 

/ 


estimated  existence  time  £ran  creation  of  particle  I  to  exposure 
1 

of  film  2. 

3 

initial  guess  at  direction  of  reference  drift 

the  angle  of  the  particle  which  is  most  clockwise /counterclockwise 
in  the  relative  angular  span. 

1 

relative  drift  angle  calculated  from  governing  equations  2  and 

1 

the  numerical  average  of  TH12,  TH23,  and  TH13 

the  relative  drift  angle  of  particle  I  with  respect  to  arbitrary 
origin  (chosen  parallel  to  film  2)  (degrees). 

the  current  reference  drift  angle  being  tested 

1  2 

relative  transverse  velocity  calculated  from  equations  2  and  3 

1  3 

relative  transverse  velocity  of  particle  I  based  on  weighted  values 
of  VETR12,  VETR23,  and  VETR13  (m/s) 

absolute  transverse  velocity  of  particle  I  (m/s) 

transverse  velocity  of  particle  reference  (m/s) 

x  component  of  drift  velocity  for  the  reference 

y 

x  component  of  relative  drift  velocity  for  particle  I 

y 

the  percentage  weight  that  an  individual  drift  velocity  calcula¬ 
tion  contributes  to  the  overall  average. 


statistical  variables  used  in  VA  vs  Vy  regression 


w  w  ro 


"PEN  JET' 


T  (I) 


-  inpact  time  of  particle  I,  taken  from  virtual  origin 


X  (I)  ) 

V  (I)  > 

2  (I)  ) 

DELT 

U  (I) 

RADIUS  (I) 
PEN  (I) 
HVOL  (I) 

V  (I,  1)  ) 

V  (I,  2)  } 

V  (I,  3)  ) 

VM  (1) 

VA  (1 ,  1)  ) 
VA  (I,  2)  J 
VA  (I ,  3)  ) 


VAM  (I) 


VP  (I 
VP  (I 
VP  (I 


,  m 
,  2)  | 
,  3)  ) 


VPM  (I) 

RGRAD  (I 
RGRAD  (I 
RGRAD  (I 


-  spacial  coordinates  for  top  of  particle  I  at  any  given  time 
(and  eventually  upon  impact)  X  =  0,  Y  =  0  denote  where  charge 
axis  intersects  target  plate;  Z  =  0  is  target  surface. 

-  time  increment  (varied  so  that  particle  travels  3  mm/timestep) 

-  penetration  velocity  of  particle  I  (calculated  only  in  the 
hypervelocity  regime  e.g.  above  2  mm/psec) 

-  crater  radius  of  particle  I 

-  penetration  of  particle  I 

-  hole  volume  created  by  particle  I 


the  y  components  of  a  vector  drawn  from  the  center  of 
z 

crater  J  (through  which  particle  I  is  instantaneously  passing) 
to  the  tip  of  particle  I 


-  magnitude  of  vector  V 


the  y  components  of  vector  VA  which  is  in  itself  the  component 
z 

of  vector  V  pointing  axially  down  crater  J  (the  crater  through 
which  particle  I  is  passing) 


-  magnitude  of  vector  VA 


the  y  components  of  vector  VP  which  is  in  itself  the  component 
2  A 

of  vector  V  pointing  perpendicular  to  the  penetration  axis  of 
crater  J  (through  whi^ch  particle  I  is  passing)  . 

Note  that  VA  *  =  V 

magnitude  of  vector  VP 


'  2) 

’  2)l 

,  3)) 


the  y  components  of  a  vector  that  is  perpendicular  to  crater  J 
z 

at  the  point  of  impact  of  particle  I,  but  in  the  coordinate 
system  of  crater  J  (e.g., the  Z’  axis  points  down  the  axis  of 
crater  J) 


RGM  (I) 


-  magnitude  of  vector  RGRAD 


GRAD  (I,  1) 
GRAD  (I,  2) 
GRAD  (I,  3) 

GM  (I) 

AR  (I,  1)  ) 
AR  (I,  2)  > 
AR  (I,  3)) 

PR  (I,  1)  1 
PR  Cl,  2)  J 
PR  (I,  3)  ) 

NCC  (I) 

econst 

ZSUBO 

RHOJ 

RHOT 

SIGJ 

SIGT 

GAMMA 

JMIN 

PENSUM 
NQ  =  0 
NQ  =  1 

NQ  »  2 

KZSUM 


the  y  components  of  a  vector  that  is  perpendicular  to  crater  J 
z 

at  th^  point  of  impact  of  particle  I,  and  in  laboratory 
coordinate  system  te.g.,  the  axis  points  perpendicular  to  the 
target  surface);  GRAD  becomes  the  axis  of  penetration  for 
crater  I  according  to  refracted  particle  wave  theory 

-  magnitude  of  vector  GRAD 


a  unit  vector  version  of  VA 


a  unit  vector  version  of  VP 


-  equals  1  when  particle  I  is  particulate;  equals  0  when  I  is 
continuous 

-  energy  constant  used  for  calculating  crater  volume  (J/mm3) 

-  standoff  (mm)  from  virtual  origin  to  target  surface 

-jet  density  P  (g/mm3) 

-  target  density  P^(g/mm  ) 

-  jet  strength  (MPa) 

-  target  strength  (MPa) 

-  (Wb 

-  the  last  crater  number  that  particle  I  has  been  known  to 
successfully  have  passed  through 

-  the  penetration  resulting  from  PENJET  calculations 

-  particle  I  is  travelling  down  hole  profile 

-  particle  I  has  just  struck  target;  t  is  decremented  until 
precise  impact  time  is  revealed 

-  particle  I  has  drifted  so  far  off  course  as  to  preclude  entry 
into  any  part  of  the  hole  profile;  it  therefore  strikes  the 
target  surface 

-  number  of  standoffs  at  which  PENJET  has  just  analyzed  a  drift 
distribut ion 


PENMAX  (K)  -  the  greatest  possible  penetration  at  standoff  #K 


an  array  of  the  standoffs  at  which  PENJET  has  just  been 
penetration  predicted  by  PENJET  at  standoff  MC 


DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


12  Commander 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22314 

2  HQDA  (DAMA-ZA;  DAMA-RA,  T.  Rorstad) 
Washington,  DC  20310 

1  Assistant  Secretary  of  the 
Army  (R5D) 

ATTN:  Assistant  for  Research 
Washington,  DC  20310 

1  Deputy  Under  Secretary  of  the 
Army  (Operations  Research) 

ATTN:  SAUS-OR,  D.  C.  Hardison 
Room  2E621,  The  Pentagon 
Washington,  DC  20310 

3  Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMD-ST 
DRCCG 
DRCDE-WM 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

5  Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRDAR-TSS 

DRDAR-LCU-D-R/T.  Stevens 
DRDAR-LCU-LT/A.  King 
DRDAR-LCU-DC/J.  Pearson 
Dover,  NJ  07801 

1  Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRDAR-TDC  (Dr.  D.  Gyorog) 
Dover,  NJ  07801 

1  Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-LEP-L 
Rock  Island,  IL  61299 

79 


1  Director 

US  Army  ARRADCOM 
Benet  Weapons  Laboratory 
ATTN:  DRDAR-LCB-TI. 

Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRDAV-E 
4300  Goodfellow  Blvd 
St.  Louis,  M0  63120 

1  Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Communications  Research 
and  Development  Command 
ATTN:  DRSEL-ATDD 
Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Fort  Monmouth,  NJ  07703 

9  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-R 
DRSMI-J 

COL  C.  M.  Matthews,  Jr. 
DRSMI-VIT 
DRCPM- HD/COL  Cass 
DRCPM-HDE-M/R.  Kasucci 
DRCPM-TO 
DRSMI-YDI. 

TOW- PM/COL  Williamson 
DRSMI  -RLA/1:.  B.  Harwell 
Redstone  Arsenal,  AL  35898 

2  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CS0-0R 
Fort  Benning,  GA  31905 


DISTRIBUTION 


No.  of 
Copies 


Organization 


3  Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  DRDTA-RWL/C.  Bradley 
DRDTA-TSL 

DRDTA-RC/Dr.  J.  Jellinek 
Warren,  MI  48090 

1  Program  Manager 

XM1  Abrams  Tank  System 
ATTN:  DRCPM-GCM-S 
Warren,  MI  48090 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 

2  Commander 

US  Army  Materials  and  Mechanics 
Research  Center 
ATTN:  DRXMR-PDD/J.  Prifti 
Tech  Lib 

Watertown,  MA  02172 

1  Commander 

US  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
NC  27709 

2  Chief  of  Naval  Research 
Department  of  the  Navy 
ATTN:  Code  427 

Code  470 

Washington,  DC  20360 

1  Commander 

Naval  Sea  Systems  Command 
Washington,  DC  20360 


Commander 

Naval  Air  Systems  Command 
ATTN:  Code  AIR-310 
Code  AIR-350 
Washington,  DC  20360 


Organization 


2  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  DG-50 

DX-21 ,  Lib  Br 

Dahlgren,  VA  22448  22448 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  730,  Lib 
Silver  Spring,  MD  20910 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  N.  Coleburn 
Silver  Spring,  MD  20910 

1  Commander 

Naval  Weapons  Center 
ATTN:  Dr.  L.  Smith 
China  Lake,  C A  93555 

2  Commander 

Naval  Weapons  Center 
ATTN:  Code  4 OS 7 

Code  3431,  Tech  Lib 
China  Lake,  CA  93555 

1  Commander 

Naval  Research  Laboratory 
Washington,  DC  20375 

1  USAF/AFRDDA 

Washington,  DC  20311 

1  AFSC/5DZ/SDOA 
Andrews  AFB 
Washington,  DC  20334 

2  US  Air  Force  Academy 
ATTN:  Code  FJS-41  (NC) 

Tech  Lib 

Colorado  Springs,  CO  80840 

1  AFATL/DLJW  (Dr.  J.  Foster) 
Eglin  AFB,  FL  32542 

1  AFWL/SUL 

Kirt land  AFB,  NM  87117 


No.  of 
Copies 


Organization 


DISTRIBUTION  LIST 

No.  of 
Copies 


Organization 


1  AFATL/DLYV  (Mr.  J.  Smith) 

Eglin  AFB,  FL  32S42 

1  AFWL/SUL  (LT  Tennant) 

Kirtland  AFB,  NM  87117 

1  AFAL 

Wright-Patterson  AFB,  OH  45433 
1  AFLC/M4WMC 

Wright-Patterson  AFB,  OH  45433 

1  Director 

Lawrence  Livermore  National  Lab. 
ATTN:  Tech  Lib 
P.  0.  Box  808 
Livermore,  CA  94550 

2  Director 

Lawrence  Livermore  National  Lab. 
ATTN:  Dr.  C.  Godfrey 
Mr.  M.  Wilkins 
P.  0.  Box  808 
Livermore,  CA  94550 

1  AFELM,  The  Rand  Corporation 
ATTN:  Library  D 
1700  Main  Street 
Santa  Monica,  CA  90406 

1  Sandia  Laboratories 
ATTN:  Tech  Lib 
Albuquerque,  NM  87115 

1  Battelle-  Memorial  Institute 
ATTN:  Mr.  Joseph  E.  Backofen 
505  King  Avenue 
Columbus,  OH  43201 

1  Dyna  East  Corporation 
ATTN:  P.  C.  Chou 
227  Hemlock  Road 
Wynnewood,  Pa  19096 


Firestone  Defense  Research 
and  Products  _ 

Division  of  the  Firestone 
Tire  and  Rubber  Company 
ATTN:  R.  Berus,  L.  Swabley 
1200  Firestone  Parkway 
Akron,  OH  44317 

Honeywell  Inc. 

Government  and  Aerospace 
Products  Division 
ATTN:  C.  R.  Hargreaves 
600  Second  Street 
Hopkins,  NM  55343 

Physics  International 
ATTN:  Ron  E.  Brown 
San  Leandro,  CA  94577 

Drexel  Institute  of  Technology 
32nd  and  Chestnut  Streets 
Philadelphia,  PA  19104 

Systems,  Science  6  Software 
ATTN;  Dr.  R.  Sedgwick 
P.  0.  Box  1620 
La  Jolla,  CA  92037 


Aberdeen  Proving  Ground 


Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 
DRXSY-G,  Mr.  Kramer 
Cdr,  USATEC0M 

ATTN:  DRSTE-TO-E 
Dir,  USACSL,  Bldg.  E3516,  EA 
ATTN:  DRDAR-CLB-PA 
DRDAR-CLN 
DRDAR-CLJ-L 


81 


USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports . 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
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